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ABSTRACT 

We discuss the X-ray properties of the cooling flows in a sample of thirty highly X- 
ray luminous clusters of galaxies observed with the ASCA and ROSAT satellites. We 
demonstrate the need for multiphase models to consistently explain the spectral and 
imaging X-ray data for the clusters. The mass deposition rates of the cooling flows, 
independently determined from the ASCA spectra and ROSAT images, exhibit good 
agreement and exceed 1000 Mq yr _1 in the largest systems. We confirm the presence 
of intrinsic X-ray absorption in the clusters using a variety of spectral models. The 
measured equivalent hydrogen column densities of absorbing material are sensitive to 
the spectral models used in the analysis, varying from average values of a few 10 20 
atom cm -2 for a simple isothermal emission model, to a few 10 21 atom cm~ 2 using our 
preferred cooling flow models, assuming in each case that the absorber lies in a uniform 
foreground screen. The true intrinsic column densities are likely to be even higher 
if the absorbing medium is distributed throughout the clusters. We summarize the 
constraints on the nature of the X-ray absorber from observations in other wavebands. 
A significant part of the X-ray absorption may be due to dust. 
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1 INTRODUCTION 

X-ray observations of clusters of galaxies show that in the 
central regions of most clusters the cooling time of the in- 
tracluster gas is significantly less than a Hubble time {e.g. 
Edge, Stewart & Fabian 1992; White, Jones & Forman 1997; 
Peres et al. 1998). This cooling is thought to lead to a slow 
net inflow of material towards the cluster centre; a process 
known as a cooling flow (see Fabian 1994 for a review). X-ray 
imaging data show that gas typically 'cools out' through- 
out the central few tens to hundreds of kpc in clusters, with 
M(r) °$ r, where M(r) is the integrated mass deposition rate 
within radius r {e.g. Thomas, Fabian & Nulsen 1987). Spa- 
tially resolved X-ray spectroscopy has confirmed the pres- 
ence of distributed cooling gas in cooling flows, with a spa- 
tial distribution and luminosity in good agreement with the 
predictions from the imaging data and cooling-flow models 
{e.g. Allen & Fabian 1997). 

For some years, the primary uncertainty with the stan- 
dard model for cooling flows was the fate of the cooled mat- 
ter {e.g. see Fabian, Nulsen & Canizares 1991). However, the 
discovery of large column densities of intrinsic X-ray absorb- 
ing material associated with cooling flows observed with the 
Solid State Spectrometer (SSS) on the Einstein Observatory 
(White et al. 1991; Johnstone et al. 1992) opened one inter- 



esting possibility. Follow-up spatially-resolved X-ray spec- 
troscopy with the Position Sensitive Proportional Counter 
(PSPC) on ROSAT {e.g. Allen et al. 1993; Irwin & Sarazin 
1995; Allen & Fabian 1997) confirmed the presence of in- 
trinsic X-ray absorbing material in cooling-flow clusters and 
further showed this material to be distributed throughout, 
but centrally-concentrated within, the cooling flows. The X- 
ray data thus identify the intrinsic X-ray absorbing material 
as a plausible sink for the cooled gas deposited by the cooling 
flows. 

In this paper we examine the X-ray properties of the 
cooling flows in a sample of 30 of the most X-ray lumi- 
nous (LboI > 10 45 erg s _1 ) clusters of galaxies known. Using 
ASCA spectra and ROSAT High Resolution Imager (HRI) 
data, we present independent determinations of the mass de- 
position rates in the cooling flows and measure the column 
densities of intrinsic X-ray absorbing material associated 
with these systems. We show that the cooling flow model 
provides a consistent description of the X-ray imaging and 
spectral data and that the observed masses of intrinsic X- 
ray absorbing material are in reasonable agreement with the 
masses expected to have been accumulated by the cooling 
flows over their lifetimes. We also discuss the evidence for 
mass deposition from the cooling flows in other wavebands. 

This is the final paper in a series that has examined 
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the impact of cooling flows on mass measurements (Allen 
1998), the fcT x - L Bo i relation (Allen & Fabian 1998a) and 
metallicity measurements (Allen & Fabian 1998b) for X-ray 
luminous clusters. In this paper, we describe the data reduc- 
tion and analysis procedures used in these works. A large 
number of other, previous studies have also analysed one or 
more of data sets included here (see references in the papers 
listed above), although these studies have not, in general, 
examined the properties of the cooling flows in the clus- 
ters. Exceptions are the deprojection analysis of Peres et al. 
(1998), and the combined ASCA/ROSAT studies of Allen 
et al. (1996), Schindler et al. (1997), Bohringer et al. (1998) 
and Rizza et al. (1998). 

The cosmological parameters Ho=50 kms -1 Mpc -1 , 
fi = 1 and A = are assumed throughout. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample selection 

Our sample was identified from clusters contained in 
the ROSAT X-ray Brightest Abell-type Cluster Sample 
(XBACS; Ebeling et al. 1996) and ROSAT Brightest Cluster 
Sample (BCS; Ebeling et al. 1998) with X-ray luminosities 
in the 0.1 — 2.4 keV band exceeding 10 45 ergs -1 . From 
this initial list, we selected for study those targets with 
ASCA X-ray spectra and ROSAT High Resolution Imager 
(HRI) data available on the Goddard Space Flight Centre 
(GSFC) public archive as of 1997 July 15. We supplemented 
this sample with a number of other, southern X-ray lumi- 
nous clusters, known to exhibit strong gravitational lens- 
ing effects (PKS0745-191, RXJ1347.5-1145, MS2137.3-2353, 
AC114) and the distant, luminous cooling-flow clusters Abell 
1068, 1704 and IRAS 09104+4109. 

The primary goal of this project was to investigate the 
X-ray properties of cooling flows in X-ray luminous clusters 
and their effects on the integrated X-ray properties of their 
host systems. We have not included either the Perseus (Abell 
426) or Coma (Abell 1656) clusters in our study, although 
these systems met our selection criteria, since both are too 
close (z = 0.0183 and 0.0232, respectively) to allow their 
integrated cluster properties to be studied using the tech- 
niques used in this paper. Detailed analyses of these clusters 
and other nearby cooling flows are discussed by Allen et al. 
(1999). 

The clusters Abell 370, 115 and 1758 also met the selec- 
tion criteria but were not included in our study since ROSAT 
HRI images show them not be single, coherent structures but 
to consist of a number of merging subunits. The integrated 
X-ray spectra for such clusters will relate to the virial prop- 
erties of the individual subclusters rather than the systems 
as a whole, and the assumptions of spherical symmetry and 
hydrostatic equilibrium required by the X-ray modeling will 
not apply. The X-ray images for the other clusters included 
in the sample do not exhibit any dramatic substructure that 
would suggest these assumptions to be invalid. Our final 
sample consists of 30 clusters, spanning the redshift range 
0.056 < z < 0.451, with a mean redshift of 0.21. 



2.2 The ASCA observations 

The ASCA (Tanaka, Inoue & Holt 1994) observations were 
made over a three-and-a-half year period between 1993 
April and 1996 December. The ASCA X-ray Telescope array 
(XRT) consists of four nested-foil telescopes, each focussed 
onto one of four detectors; two X-ray CCD cameras, the 
Solid-state Imaging Spectrometers (SIS; SO and SI), and 
two Gas scintillation Imaging Spectrometers (CIS; G2 and 
G3). The XRT provides a spatial resolution of ~ 3 arcmin 
Half Power Diameter (HPD) in the energy range 0.3 — 12 
keV. The SIS detectors provide good spectral resolution 
[AE/E = 0.02(£/5.9keV)~ ' 5 ] over a 22 x 22 arcmin 2 field 
of view. The GIS detectors provide poorer energy resolution 
[AE/E = 0.08(£/5.9keV)~ - 5 ] but cover a larger circular 
field of view of ~ 50 arcmin diameter. 

For our analysis we have used the screened event lists 
from the revl processing of the data sets available on the 
GSFC ASCA archive (for a detailed description of the revl 
processing see the GSFC ASCA Data Reduction Guide, 
published by GSFC.) The ASCA data were reduced us- 
ing the FTOOLS software (version 3.6) issued by GSFC, 
from within the XSELECT environment (version 1.3). Fur- 
ther data-cleaning procedures as recommended in the ASCA 
Data Reduction Guide, including appropriate grade selec- 
tion, gain corrections and manual screening based on the 
individual instrument light curves, were followed. A sum- 
mary of the ASCA observations, including the individual 
exposure times after all screening procedures were carried 
out, is provided in Table 1. 

Spectra were extracted from all four ASCA detectors 
(except in those few cases where the SI data were lost due 
to saturation problems caused by flickering pixels in the 
CCDs). The spectra were extracted from circular regions, 
centred on the peaks of the X-ray emission. For the SIS 
data, the radii of the regions used were selected to minimize 
the number of chip boundaries crossed (thereby minimizing 
the systematic uncertainties introduced by such crossings) 
whilst covering as large a region of the clusters as possible. 
Data from the regions between the chips were masked out 
and excluded. The final extraction radii for the SIS data are 
summarized in Table 2. Also included in that Table are the 
chip modes used in the observations (whether 1,2 or 4 chip 
mode) and the number of chips from which the extracted 
data were drawn. For the GIS data a constant extraction 
radius of 6 arcmin was used. [We note that for Abell 2142, 
the 2 arcmin (3 arcmin) radius region surrounding the X- 
ray bright Seyfert-1 galaxy 1556+274 (offset by ~ 4 arcmin 
from the X-ray centre of the cluster) was masked out and 
excluded from the analysis of the SIS (GIS) data.] 

For the GIS observations, and SIS observations of clus- 
ters in regions of low Galactic column density (Ah < 5 x 10 20 
atom cm -2 ), background subtraction was carried out using 
the 'blank sky' observations of high Galactic latitude fields 
compiled during the performance verification stage of the 
ASCA mission. For such data sets, the blank-sky observa- 
tions provide a reasonable representation of the cosmic and 
instrumental backgrounds in the detectors. The background 
data were screened and grade selected in the same manner 
as the target observations and background spectra were ex- 
tracted from the same regions of the detectors as the cluster 
spectra. For the SIS observations of clusters in directions 
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of higher Galactic column density, background spectra were 
extracted from regions of the chips that were relatively free 
from foreground cluster emission. 

For the SIS data, response matrices were generated us- 
ing the FTOOLS SISRMG software. Where the spectra cov- 
ered more than one chip, response matrices were created 
for each chip, which were then combined to form a counts- 
weighted mean matrix. For the GIS analysis, the response 
matrices issued by GSFC on 1995 March 6 were used. 

2.3 The ROSAT observations 

The ROSAT HRI observations were carried out between 
1991 November and 1995 June. The HRI provides a ~ 5 
arcsec (FWHM) X-ray imaging facility covering a ~ 40 x 40 
arcmin 2 field of view (David et al. 1996). Reduction of the 
data was carried out with the Starlink ASTERIX package. 
X-ray images were created on a 2 x 2 arcsec 2 pixel scale, 
from which centres for the cluster X-ray emission were de- 
termined. Where more than one observation of a source was 
made, a mosaic was constructed from the individual obser- 
vations. For the cooling flow and intermediate clusters, the 
X-ray centres were identified from the peaks of the X-ray 
surface brightness distributions (which are easily determined 
from the HRI images). For the non-cooling flow clusters, the 
X-ray emission is not as sharply-peaked and for these sys- 
tems we have identified the X-ray centres with the results 
from iterative determinations of the centroids within a circu- 
lar aperture of fixed radius. For Abell 665, 2163 and AC114 
a 1 arcminute radius aperture was used. For Abell 2744, 
773, 2218 and 2219 a 2 arcmin aperture was adopted, and 
for Abell 520 a 3 arcmin aperture was used. For two of the 
clusters included in the sample, Abell 665 and 1413, ROSAT 
HRI images were not available at the time of writing and for 
these clusters PSPC imaging data have been used instead. A 
summary of the ROSAT observations and the X-ray centers 
for the clusters is given in Table 3. 

2.4 Classification of clusters as cooling flow (CF) 
and non-cooling flow (NCF) systems 

For the purposes of this paper, we have classified the clus- 
ters into subsamples of cooling-flow (hereafter CF) and non- 
cooling flow (NCF) systems. CFs are those clusters for which 
the upper (90 percentile) limit to the central cooling time, 
as determined from the deprojection analysis of the ROSAT 
HRI X-ray images (Section 4), is less than 10 10 yr. (The 
'central' cooling time is the mean cooling time of the cluster 
gas in the innermost bin included in the deprojection analy- 
sis, which is of variable size. The use of a fixed physical size 
of 100 kpc for the central bin leads to similar results; Sec- 
tion 4.) Using this simple classification scheme we identify 
21 CFs and 9 NCFs in our sample. The mean redshift for 
the subsamples of both CF and NCF clusters is z = 0.21. 
The fraction of CFs in our sample is 70 per cent, the same 
as that determined by Peres et al. (1998) from a study of 
a larger, flux-limited sample of clusters (primarily at lower 
redshift; z — 0.056) using a similar classification scheme. 



3 SPECTRAL ANALYSIS OF THE ASCA DATA 



3.1 The spectral models 

The modeling of the X-ray spectra has been carried out using 
the XSPEC spectral fitting package (version 9.0; Arnaud 
1996). For the SIS data, only counts in pulse height analyser 
(PHA) channels corresponding to energies between 0.6 and 
10.0 keV were included in the analysis (the energy range 
over which the calibration of the SIS instruments is best- 
understood). For the GIS data only counts in the energy 
range 1.0 — 10.0 keV were used. The spectra were grouped 
before fitting to ensure a minimum of 20 counts per PHA 
channel, allowing x 2 statistics to be used. 

The spectra have been modeled using the plasma codes 
of Kaastra &: Mewe (1993; incorporating the Fe L calcula- 
tions of Liedahl, Osterheld & Goldstein 1995, in XSPEC 
version 9.0) and the photoelectric absorption models of 
Balucinska-Church & McCammon (1992). The data from all 
four ASCA detectors were analysed simultaneously (except 
where the SI data were lost due to chip saturation problems) 
with the fit parameters linked to take the same values across 
the data sets. The exceptions to this were the emission mea- 
sures of the ambient cluster gas in the four detectors which, 
due to the different extraction radii used and residual un- 
certainties in the flux calibration of the instruments, were 
allowed to fit independently. 

The spectra were examined with a series of spectral 
models. Model A, consisted of an isothermal plasma in col- 
lisional equilibrium, at the optically-determined redshift for 
the cluster, and absorbed by the nominal Galactic column 
density (Dickey & Lockman 1990). The free parameters in 
this model were the temperature (kT) and metallicity (Z) 
of the plasma and the emission measures in the four detec- 
tors. (The metallicities are measured relative to the solar 
values of Anders & Grevesse (1989), with the different el- 
ements assumed to be present in their solar ratios.) The 
second model, model B, was the same as model A but with 
the absorbing column density (A/h) also included as a free 
parameter in the fits. The third model, model C, was the 
same as model A but with an additional component explic- 
itly accounting for the emission from the cooling flows. The 
material in the cooling flows is assumed to cool at constant 
pressure from the ambient cluster temperature, following the 
prescription of Johnstone et al. (1992). The normalization 
of the cooling-flow component was parameterized in terms 
of a mass deposition rate, Ms, which was a free parameter 
in the fits. The mass deposition rate is linked to take the 
same value in all four detectors, scaled by a normalization 
factor proportional to the total flux measured in that detec- 
tor. The metallicity of the cooling gas was assumed to be 
equal to that of the ambient ICM. The emission from the 
cooling flows was also assumed to be absorbed by an intrin- 
sic column density, A Ah, of cold gas which was a further 
free parameter in the fits. The abundances of metals in the 
absorbing material were fixed to their solar values (Anders 
& Grevesse 1989) although the effects of varying these val- 
ues were examined. Finally, a fourth spectral model, model 
D, was examined which was similar to model C but with 
the constant-pressure cooling flow replaced with a second 
isothermal emission component. The temperature and nor- 
malization of this second emission component were included 
as free parameters in the fits. (As with model C, the nor- 
malizations were linked to the same value in all four detec- 
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tors, scaled by the appropriate normalization factors.) The 
second emission component was also assumed to be intrinsi- 
cally absorbed by a column density of cold gas, ANh, which 
was a free parameter in the fits. The metallicities of the two 
emission components were linked to take the same values. 

Fig. 1 shows the ASCA data and best-fitting spectral 
models for four of the clusters in the sample; Abell 1704, 
2029, 2204 and 2219. Table 4 summarizes the fit results for 
all of the clusters using the four spectral models. For each 
cluster we list the best-fitting parameter values and 90 per 
cent (Ax 2 = 2.71) confidence limits. The mass deposition 
rates (Ms), 2—10 keV X-ray luminosities (Zoc) and bolo- 
metric luminosities (Lboi) are the values measured in the 
G3 detector. Note that for the clusters at relatively low red- 
shifts (z < 0.1), we have corrected the Lx and L Bo i values 
for the emission arising at radii > 6 arcmin by scaling from 
the David et al. (1993) results, which were based on GINGA 
and Einstein MPC observations. 

3.2 The goodness-of-fit 

Table 5 summarizes the goodness-of-fit measurements ob- 
tained with the different spectral models. The tabulated re- 
sults are the probabilities of exceeding the \ 2 values ob- 
tained, assuming in each case that the model correctly de- 
scribes the spectral properties of the clusters. Goodness-of- 
fit values < 10~ 2 may be regarded as indicating a formally 
unacceptable fit to the data. 

The results from the spectral analysis show that model 
A (the isothermal model with the absorbing column density 
fixed at the nominal Galactic value) provides a reasonable 
description of most (21/30) of the clusters. However, model 
A fails significantly for the massive CF clusters Abell 1068, 
1795, 2029 and 2204. The significance of the improvements 
to the fits obtained with spectral model B over model A 
have been evaluated using the F-test for the introduction of 
an additional free parameter (Bevington 1969). These sig- 
nificances are listed in Column 6 of Table 5. We see that the 
improvements obtained when including the column density 
as a free parameter in the fits with the isothermal models 
are significant at > 90 per cent confidence for 23 of the 30 
clusters (and would also be highly significant for Abell 478 if 
the nominal Dickey & Lockman 1990 value for the Galactic 
column density, rather than the X-ray value determined by 
Allen & Fabian 1997 were used in model A). Model B pro- 
vides a statistically acceptable fit for 27 of the 30 clusters. 

The results on the goodness-of-fit determined with spec- 
tral models C and D are summarized in columns 4 and 
5 of Table 5. The only clusters that these models do not 
adequately describe are the two nearest, brightest systems; 
Abell 1795 and Abell 2029. Spectral models C and D include 
an absorption component acting on the entire cluster emis- 
sion which is normalized to the appropriate Galactic column 
densities determined from HI studies (Dickey & Lockman 
1990). We have examined the significance of the improve- 
ments to the fits obtained when allowing the Galactic col- 
umn density to also be a free parameter in the fits with 
these models. The fits were not significantly improved for 
any of the clusters except PKS0745-191, for which the Galac- 
tic column density is most uncertain. (The Galactic column 
density to PKS0745-191 varies between 3.2 and 4.4 x 10 21 
atom cm~ 2 in the surrounding regions studied by Dickey 



& Lockman 1990). We find A\ 2 = 7.9 with model C and 
Ax 2 = 20.2 with model D for PKS0745-191 (indicating im- 
provements significant at > 99 per cent confidence) when 
allowing the Galactic column density to be free. The pre- 
ferred value for the Galactic column density for this cluster 
from the ASCA data is ~ 3.5 x 10 21 atom cm~ 2 . For Abell 
1795, we also find a marginal improvement to the fits with 
models C and D when allowing the Galactic column den- 
sity to be a free parameter (A^ 2 = 4.5 with Model C and 
Ax = 5.4 for model D, which is significant at the 96.5 
and 98 per cent level, respectively) with a preferred value of 
< 8 x 10 19 atom cm~ 2 . The agreement between the Galactic 
column densities inferred from the 21 cm and X-ray data for 
most of the clusters with spectral models C and D provides 
further support for the validity of these models. 

3.3 The requirement for multiphase models 

The statistical improvements obtained by introducing an 
additional emission component into the fits with the single- 
temperature model (model B) are summarized in Columns 7 
and 8 of Table 5. We list the results obtained for both cases, 
where the extra emission component is modeled as a cooling 
flow (for the CF systems only) or a second isothermal emis- 
sion component. (Essentially we use spectral models C and 
D, but also include the Galactic column density as a free 
parameter it the fits so as to permit a direct comparison 
with the goodness of fit measured with spectral model B). 
The significance of the improvements to the fits obtained 
with the multiphase models with respect to model B have 
been quantified using an F-test for the introduction of 2 
or 3 extra parameters (for models C and D respectively). 
We find that the improvements with model D are signifi- 
cant at > 90 per cent confidence for 21 of the 30 clusters. 
The only clusters for which we find no significant improve- 
ment with the introduction of a second emission compo- 
nent are the NCF clusters Abell 520, 773, 2218 and AC114, 
and the more distant and/or less- luminous CF clusters Abell 
963, MS1358.4+6245, MS1455.0+2232, MS2137.3-2353 and 
Abell 2390, which have amongst the lowest numbers of 
counts in their ASCA spectra. For the CF clusters (the only 
clusters to which model C was applied) the improvements 
obtained with the cooling-flow model over model B are sig- 
nificant at > 90 per cent confidence in 13 of 21 systems, with 
again those CF clusters with the lowest numbers of total 
counts in their ASCA spectra exhibiting the least-significant 
improvements. 

We thus find that most CF clusters, and a number of 
NCF systems (in general those systems with the best signal- 
to-noise ratios in their ASCA spectra) exhibit significant 
improvements to their fits with the use of multiphase, over 
single-phase, models. These improvements generally indicate 
the presence of emission from material cooler than the mean 
cluster temperatures (see Table 4). For the CF clusters, this 
can be naturally understood as being due to emission as- 
sociated with the cooling flows. For the NCF clusters, the 
improvements to the fits obtained with model D are likely to 
reflect the presence of merging subclusters with lower virial 
temperatures. This is supported by optical, radio and X-ray 
imaging (e.g. Edge et al. 1992; Buote & Tsai 1996; Marke- 
vitch 1996; Feretti, Giovannini & Bohringer 1997; Rizza et 
al. 1998) and gravitational lensing studies (e.g. Kneib et al. 
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Figure 1. (Upper panels) The ASCA data and best-fitting models and (lower panels) residuals to the fits (in units of x) for a representative 
subsamplc of the clusters included in the study. Clockwise, from top left, we show Abell 2029, 2219, 1704 and 2204. For the 3 CF clusters 
(Abell 2029, 2204 and 1704), the results determined with spectral model C are shown. For the NCF cluster, Abell 2219, spectral model 
B is shown. For clarity, only the SO and G3 data are plotted. The figure illustrates the range in quality of the ASCA data included in 
the study. 



1995; Smail et al. 1995, 1997; Squires et al. 1997; Allen 1998) 
of the NCF clusters, which show that these systems typically 
exhibit complex morphologies and centroid shifts indicative 
of merger events. 



3.4 Two-temperature versus cooling-flow models 

The results in Table 4 show that spectral model D typically 
provides at least as good a fit to the ASCA spectra for the 
CF clusters as the cooling-flow model (model C). This has 
sometimes been taken to indicate that the X-ray gas in these 



systems is distinctly two-phase, with the cooler phase being 
due to the dominant cluster galaxy (e.g. Makishima 1997; 
Ikebe et al. 1999). However, when interpreting these results 
it is important to recall that the two-temperature model 
provides a more flexible fitting parameterization, with an 
extra degree of freedom, when applied to ASCA observa- 
tions. Simulated cooling-flow spectra constructed with spec- 
tral model C, including plausible levels of intrinsic absorp- 
tion (Section 6) and observed at the spectral resolution and 
signal-to-noise levels typical of ASCA observations, are in- 
variably well-described by two-temperature models. This is 
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Energy (keV) 

Figure 2. (Upper panel) Simulated ASCA SIS spectrum of a 7 
keV cluster with a metallicity of OAZq and a cooling flow, in- 
trinsically absorbed by a column density of 4 X 10 21 atom cm -2 , 
accounting for 30 per cent of the total 2 — 10 keV flux. A count 
rate of let s _1 , an exposure time of 40ks and a Galactic col- 
umn density of 10 20 atom cm -2 have been assumed. The best- 
fitting two-temperature model (which provides a \ 2 of 263 for 
245 degrees of freedom) is overlaid. The best-fit parameters are 
fcTi = 6.7 ± 0.4 keV, kT 2 = 1.4 ± 0.2 keV, Z = 0.41 ± O.O4Z , 
with an intrinsic column density acting on the cooler component 
of 7.6 ± 1.8 x 10 21 atom cm ). (Lower panel) The residuals to 
the fit in units of count s -1 keV _1 . 

illustrated in Fig. 2 where we show a simulated ASCA SIS 
spectrum for a 7 keV cluster with a metallicity of OAZq, 
containing a cooling flow (intrinsically absorbed by a col- 
umn density of 4 x 10 21 atom cm~ 2 ) accounting for 30 per 
cent of the total 2 — 10 keV flux. A count rate of let s _1 , 
an exposure time of 40ks and a Galactic column density of 
10 20 atom cm have been assumed. Overlaid, we show the 
best-fitting two-temperature model, which provides a good 
fit to the simulated data. 

Thus, even in the case where the constant-pressure cool- 
ing flow model provides an exact description of the data 
(i.e. in the simulations), the two-temperature model pro- 
vides a similarly good fit. In a real cluster, where spectral 
model C will undoubtedly over-simplify the true situation, 
the more flexible two-temperature model (model D) is likely 
to provide a better match to the observations. Thus, our 
finding that spectral model D often provides a better fit to 
the observed cluster spectra than model C, only shows that 
the constant pressure cooling flow model over-simplifies the 
true spectra of the cooling flows in the clusters (see also Sec- 
tion 6.4). We note that such tendencies would be enhanced 
by the presence of additional cool components in the spec- 
tra, for example due to temperature gradients at large radii 
(e.g. Markevitch et al. 1998). Finally, we note that the very 
high luminosities associated with the cool emission compo- 
nents in clusters like RXJ1347. 5-1145 (with L coo i ~ a few 
10 45 ergs _1 ) are difficult to explain with models in which 



the cool emission is due to the interstellar medium of the 
central cluster galaxies. 

4 DEPROJECTION ANALYSIS OF THE 
X-RAY IMAGES 

4.1 Method and primary results 

The analysis of the HRI imaging data has been carried out 
using an extensively updated version of the deprojection 
code of Fabian et al. (1981). Azimuthally-averaged X-ray 
surface brightness profiles were determined for each clus- 
ter from the HRI images. These profiles were background- 
subtracted, corrected for telescope vignetting and re-binned 
to provide sufficient counts in each annulus to allow the anal- 
ysis to be extended to radii of at least 500 kpc. (Bin sizes of 
8-24 arcsec were used.) 

With the X-ray surface brightness profiles as the pri- 
mary input, and under assumptions of spherical symmetry 
and hydrostatic equilibrium, the deprojection technique can 
be used to study the basic properties of the intracluster gas 
(temperature, density, pressure, cooling rate) as a function 
of radius. The deprojection code uses a monte-carlo method 
to determine the statistical uncertainties on the results and 
incorporates the appropriate HRI spectral response matrix 
issued by GSFC. The cluster metallicities were fixed at the 
values determined from the spectral analysis (Table 4). The 
absorbing column densities were fixed at the appropriate 
Galactic values (Table 1). 

The deprojection code requires the total mass profiles 
for the clusters (which define the pressure profiles) to be 
specified. We have iteratively determined the mass profiles 
that result in deprojected temperature profiles (which ap- 
proximate the mass-weighted temperature profiles in the 
clusters) that are isothermal within the regions probed by 
the HRI data (the central 0.5 — 1 Mpc) and which are con- 
sistent with the spectrally-determined temperatures (Sec- 
tion 3). The assumption of approximately isothermal mass- 
weighted temperature profiles in the central regions of the 
clusters is supported by the following evidence: firstly, ASCA 
observations of nearby cooling flows show that in the cen- 
tral regions of these systems the gas is multiphase, but that 
the bulk of the X-ray gas there has a temperature close to 
the cluster mean (e.g. Fabian et al. 1994b; Fukazawa et al. 
1994; Matsumoto et al. 1996; Ikebe et al. 1999). Secondly, 
combined X-ray and gravitational lensing studies of CF clus- 
ters (e.g. Allen 1998) show that approximately isothermal 
mass- weighted temperature profiles in the cores of such clus- 
ters lead to excellent agreement between their X-ray and 
gravitational-lensing masses. Thirdly, the use of approxi- 
mately constant mass-weighted temperature profiles implies 
a more plausible range of initial density inhomogencities in 
the clusters than would be the case if the temperature pro- 
files decreased within their cores (Thomas, Fabian & Nulsen 
1987). Finally, the use of approximately isothermal mass- 
weighted temperature profiles in the deprojection analyses 
leads to independent determinations of the mass deposition 
profiles in the cooling flows, from the X-ray spectra and 
imaging data, in excellent agreement with each other (e.g. 
Allen & Fabian 1997). We note that the assumption of a con- 
stant mass- weighted deprojected temperature profile is con- 
sistent with measurements of a decreasing emission-weighted 
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temperatures in the cores of many CF clusters (e.g. Waxman 
& Miralda-Escude 1995). 

The mass profiles for the clusters were parameterized as 
isothermal spheres (Equation 4-125 of Binney & Tremaine 
1987) with adjustable core radii (r c ) and velocity disper- 
sions (a). The core radii were adjusted until the tempera- 
ture profiles determined from the deprojection code became 
isothermal. The velocity dispersions were then adjusted un- 
til the temperatures determined from the deprojection code 
came into agreement with the spectrally-determined values. 
Errors on the velocity dispersions are the range of values 
that result in isothermal deprojected temperature profiles 
that are consistent, at the 90 per cent confidence limit, with 
the spectrally-determined temperature results. Estimates of 
the thermal gas pressure in the outermost radial bins used 
in the analysis are also required by the deprojection code 
and were determined iteratively. (The uncertainties on the 
outer pressure estimates do not significantly affect the re- 
sults presented here.) Although the deprojection method of 
Fabian et al. (1981) is essentially a single-phase technique, 
we note that it produces results in good agreement with 
more detailed multiphase treatments (Thomas, Fabian & 
Nulsen 1987) and, due to its simple applicability at large 
radii in clusters, is better-suited to the present project. 

The mass distributions determined from the deprojec- 
tion analysis are summarized in columns 3 and 4 of Table 
6. These simple parameterizations permit direct compar- 
isons with independent mass constraints from gravitational 
lensing and dynamical studies (e.g. Allen 1998). We note 
that for a few of the clusters (primarily the brightest CFs 
with the best data) the single-component mass models can- 
not adequately satisfy the requirement for approximately 
isothermal deprojected temperatures profiles and for these 
systems a significantly better match was obtained by intro- 
ducing a second 'linear' mass component, truncated at a 
specified outer radius. The clusters requiring the additional 
mass components were Abell 586 (2 x 10 11 Mokpc -1 within 
the central 20 kpc), Abell 2029 (3 x 10 11 M Q kpc _1 within 
the central 20 kpc), Abell 2219 (10 11 Mokpc" 1 within the 
central 30 kpc), Abell 478 (4 x 10 11 M Q kpc _1 within the 
central 20 kpc and an adjustment of the cluster parameters 
to a = 840igo kms" 1 and r c = 100 kpc) and Abell 1795 
(10 11 Mokpc -1 within the central 20 kpc, and an adjust- 
ment of the cluster parameters to a = 720^20 km s _1 and 
r c = 60 kpc). 

The core radii of the mass distributions determined 
from the deprojection analysis are similar to the values mea- 
sured from simple '/?— model' fits to the X-ray surface bright- 
ness profiles (e.g. Cavaliere & Fusco-Femiano 1976, 1978). 
It is well known that the /3-model does not provide a good 
match to the X-ray surface brightness profiles of the bright- 
est CF clusters when their central regions are included in 
the analysis (e.g. Jones & Forman 1984). However, where 
/3-models are used to estimate the mass core radii in CF 
clusters, the central regions of these clusters should not sim- 
ply be excised from the fits. The central surface brightness 
peaks in CF clusters trace the presence of sharp central den- 
sity rises in these systems and are not only the result of 
excess emission due to cooling gas. (It is the high central 
densities that lead to the high cooling rates.) Although fully 
accounting for the effects of cooling on the central density 
measurements is difficult with current data, simply excising 
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Figure 3. The mass deposition rates determined from the ASCA 
spectra (Mg) versus the values measured from the deprojection 
analysis of the ROSAT HRI imaging data (Mi). The dashed line 
is the best-fitting power-law model. No correction for the effects 
of intrinsic absorption on the deprojection results has been made. 



the inner regions of the clusters can lead to overestimates 
of the mass core radii and corresponding underestimates of 
the central gas and total mass densities. 

The basic results on the cooling flows from the depro- 
jection analysis are summarized in columns 5 — 8 of Table 6. 
For those clusters in common, the results obtained are gen- 
erally in good agreement with values reported from previous 
works (e.g. Edge et al. 1992, White et al. 1997; Peres et al. 
1998). 



5 THE MASS DEPOSITION RATES FROM 
THE COOLING FLOWS 

5.1 Comparison of the spectral and imaging 
results 

The spectral and image deprojection analyses discussed in 
Sections 3 and 4 provide essentially independent estimates 
of the mass deposition rates in the clusters. A comparison of 
the results obtained from these analyses therefore provides 
a test of the validity of the cooling flow model. 

The deprojection method describes the X-ray emission 
from a cluster as arising from a series of concentric spherical 
shells. The luminosity in a particular shell, j, may be written 
as the sum of four components (Arnaud 1988). 



Lj = AMjHj + AM,A$, + 



j-i 

E 

.i=l 



AMi(A$>j +AHj 



(1) 



where AM,- is the mass deposited in shell j, Hj is the en- 
thalpy of the gas in shell j, and A$j is the gravitational en- 
ergy released in crossing that shell. AM is the mass 
flow rate through shell j, and AHj the change of enthalpy 
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Figure 4. The mass deposition rates determined from the ASCA 
spectra (Ms) versus the values measured from the deprojection 
analysis, corrected for the effects of intrinsic absorption due to 
cold gas (Mc). The dashed line indicates equality between the 
values. 



of the gas as it moves through that shell. The first term 
in equation [j] thus accounts for the enthalpy of the gas de- 
posited in shell j. The second term is the gravitational work 
done on the gas deposited in shell j. The third and fourth 
terms respectively account for the gravitational work done 
on material flowing through shell j to interior radii, and the 
enthalpy released by that material as it passes through the 
shell. 

In any particular shell, the densest material in the cool- 
ing flow is assumed to cool out and be deposited. Since the 
cooling time of this material will be short compared to the 
flow time, the cooling can be assumed to take place at a 
fixed radius. Thus, the luminosity contributed by the first 
term in equation 1 should have a spectrum that can be ap- 
proximated by gas cooling at constant pressure from the 
ambient cluster temperature i.e. the same spectrum as the 
cooling component incorporated into the spectral analysis 
with model C (Section 3). 

For the bulk of the material continuing to flow inwards 
towards the cluster centre, the cooling via X-ray emission 
is assumed to be offset by the gravitational work done on 
the gas as it moves inwards. The emission accounted for in 
the second and third terms of equation 1 should therefore 
have a spectrum that can be approximated by an isothermal 
plasma at the appropriate ambient temperature for the clus- 
ter i.e. the spectrum of the isothermal emission component 
in model C. Since the mass-weighted temperature profiles in 
the clusters are assumed to remain approximately isother- 
mal with radius, the luminosity contributed by the fourth 
term of equation |l| should be negligible. 

The mass deposition rates, Mi, listed in Table 6 are the 
mass flow rates (X^-i AM<) determined from the depro- 
jection analysis, at the point where the mean cooling time 



of the cluster gas first exceeds the Hubble time (1.3 x 10 10 
yr). Thus, if the cooling flow model is correct, and the cool- 
ing flows have been undisturbed for a significant fraction of 
a Hubble time, the mass deposition rates determined from 
the deprojection analysis should be similar to those mea- 
sured independently from the spectral data. 

Fig. 3 shows the mass deposition rates determined from 
the spectral analysis (Ms) versus the image deprojection 
results (Mi), for those clusters for which both quantities 
have been measured. We see that the results exhibit an 
approximately linear correlation; a fit to the data with a 
power- law model of the form Ms = PM® , using the Akri- 
tas & Bershady (1996) bisector modification of the ordinary 
least-squares statistic, gives a best-fitting power-law slope of 
1.0 ± 0.2 (where the error is the standard deviation deter- 
mined by bootstrap re-sampling), although the spectrally- 
determined mass deposition rates typically exceed the de- 
projection results by a factor of 2 — 3. However, the spectral 
analysis presented in Section 3 also requires that the cooling 
gas is intrinsically absorbed by equivalent hydrogen column 
densities of, typically, a few 10 21 atom cm" 2 (model C). 
The mass deposition rates determined from the deprojec- 
tion analysis must therefore also be corrected for the effects 
of this absorbing material. 

The corrections for the effects of intrinsic absorption on 
the deprojection results have been carried out by re-running 
the deprojection analysis with the absorbing column densi- 
ties set to the total values determined with spectral model 
C. We assume that the absorption is due to cold gas with 
solar metallicity. If, however, the intrinsic absorption were 
instead due to dust, a possibility examined in more detail in 
Section 7, then the required correction factors would likely 
be reduced by a few tens of per cent. (The introduction 
of a simple OIK absorption edge at E ~ 0.54keV, such as 
might be associated with oxygen-rich, silicate dust grains, 
generally provides as good a description of the intrinsic ab- 
sorption detected in the 0.6 — lO.OkeV ASCA spectra as a 
cold, gaseous absorber.) Since the intrinsic column densities 
measured with model C are redshifted quantities, we set 
the total column densities used in the revised deprojection 
analysis, which assumes zero redshift for the absorber, to be 
Ah + AAh/(1 + z) 3 , where A^h is the Galactic column den- 
sity). The absorption-corrected mass deposition rates (Mc) 
so determined are summarized in Table 7. Fig. 4 compares 
the spectrally-determined mass deposition rates with the 
absorption-corrected deprojected values. We see that once 
the presence of intrinsic absorption has been accounted for 
in a consistent manner in the imaging and spectral analyses, 
the results on the mass deposition rates in the clusters show 
good agreement. 

5.2 The fraction of flux from the cooling flows 

Table 8 lists the fractions of the total 2—10 keV X-ray 
fluxes from the CF clusters associated with cooling g 
determined with spectral model C. The results range from ~ 
10 per cent for systems like Abell 1689 and 2142 to > 30 per 
cent for Zwicky 3146, Abell 1068, RXJ1347-1145 and Abell 
2204. Such a range of measurements is consistent with the 
results of Peres et al. (1998), from a deprojection analysis 
of a flux-limited sample of the brightest clusters known. 
The fraction of the X-ray flux from a cluster due to it's 
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cooling flow can be expected to increase with the age of the 
flow. It is therefore likely that the cooling flows in Abell 1689 
and 2142 are younger than those in Abell 2204 and Zwicky 
3146. We note that the X-ray data are not, in general, of 
sufficient quality to place firm constraints on the ages of the 
cooling flows in these clusters (for a discussion of cooling 
flow age measurements see Allen et al. 1999). 

It is interesting that the clusters with the largest 
cooling-flow flux fractions (with the exception of Abell 1413 
and 2261) are also amongst the most optically line-luminous 
cooling flows known (e.g. Crawford et al. 1999). The pro- 
duction of powerful optical emission lines and associated 
UV/blue continuum emission is known to be due, at least in 
part, to the formation of massive stars at the centres of cool- 
ing flows (e.g. Johnstone, Fabian & Nulsen 1987; Allen 1995; 
Cardiel et al. 1995, 1998; McNamara et al. 1996; Voit & Don- 
ahue 1997). The results presented here are consistent with 
the idea that such star formation is typical of undisturbed 
cooling flows. It is also interesting to note that the clusters 
with the largest cooling flow flux fractions and most opti- 
cally line-luminous emission line nebulae also tend to have 
the shortest cooling times averaged over their central 100 
kpc regions, with £ioo S; 3Gyr (see also Peres et al. 1998). 



6 INTRINSIC X-RAY ABSORPTION IN 
CLUSTERS 

6.1 The measured column densities and model 
dependency of the results 

The equivalent column densities of intrinsic X-ray absorbing 
material inferred from the spectral analysis are summarized 
in Table 9. All three of the spectral models incorporating a 
free-fitting absorption component (models B,C,D) reveal the 
presence of excess X-ray absorption over and above the nom- 
inal Galactic values along the lines of sight to the clusters 
(Dickey & Lockman 1990), although the measured column 
densities are sensitive to the spectral model used. 

Using the simple isothermal model with free-fitting ab- 
sorption (model B), we determine a mean excess column 
density for the whole sample of 30 clusters of 3.8 ±3.2 x 10 20 
atom cm~ 2 . When we examine the CF and NCF clusters 
separately, we find that both subsamples exhibit intrinsic 
absorption, at a similar level on average. For the CF clus- 
ters, the mean excess column density measured with spec- 
tral model B is 3.5 ± 3.0 x 10 20 atom cm" 2 . For the NCF 
systems, the value is 4.4 ± 3.9 x 10 20 atom cm -2 . The ap- 
plication of a Students t-test (accounting for the possibility 
of unequal variances in the two distributions; Press et al. 
1992) shows the mean excess column densities for the CF 
and NCF clusters, determined with spectral model B, to dif- 
fer at only the 45 per cent confidence level. The application 
of a Kolmogorov-Smirnov test shows the two subsamples to 
be drawn from different parent populations at only the ~ 10 
per cent confidence level. Thus, our analysis with spectral 
model B suggests that intrinsic absorption is not confined 
to CF clusters. 

The results on the column densities are quite different, 
however, when the more-sophisticated multiphase spectral 
models are used. As discussed in Section 3, the multiphase 
models generally provide a better description of the X-ray 



properties of the CF clusters (and some of the brighter NCF 
systems). Using spectral model C, our preferred model for 
the CF clusters in that it provides a consistent description of 
the spectral and imaging X-ray data for these systems (Sec- 
tion 5.1), we determine a mean intrinsic column density act- 
ing on the cooling-flow components of 3.4 ±1.3 x 10 21 atom 
cm~ 2 . Using the two-temperature spectral model (model 
D), we determine a mean intrinsic column density acting 
on the cooler emission components in the CF systems of 
6.3 ±4.7 x 10 21 atom cm~ 2 . Thus, the intrinsic column den- 
sities for the CF clusters measured with the multiphase mod- 
els (models C and D) are similar and approximately an or- 
der of magnitude larger than the values inferred using the 
single-phase model B. These results demonstrate the need 
for adopting an appropriate spectral when attempting to 
measure the column densities of absorbing material in CF 
clusters. We also note that where spectral model D provided 
a significant improvement to the fits to the NCF clusters, 
with respect to model B, the measured column densities 
acting on the cooler emission components were also signif- 
icantly larger than the values determined using the single- 
phase models. 

We have also measured the intrinsic column densities 
in the CF clusters using one further spectral model, referred 
to in Table 9 as model C. Model C is identical to model 
C except that it assumes that the excess absorption acts 
on the entire cluster spectrum, rather than just the cooling 
gas, and that the absorbing material lies at zero redshift. 
Model C has been used in a number of previous studies and 
is included here for comparison purposes. The mean excess 
column density for the CF clusters determined with spectral 
model C is 7.2 ± 3.9 x 10 20 atom cm -2 . 

Allen & Fabian (1997) present results from an X-ray 
colour deprojection study of 18 clusters observed with the 
ROSAT PSPC. These authors determine intrinsic column 
densities across the central 30 arcsec (radius) regions of 
Abell 478, 1795 and 2029, using spectral model C, of 
7.40 ± 0.63, 2.20 ± 0.20 and 0.83 ± 0.42 x 10 20 atom cm" 2 , 
respectively. The results for Abell 478 and Abell 1795 are 
in excellent agreement with those presented here, although 
the Allen & Fabian (1997) value for Abell 2029 is ~ 4 times 
smaller than our result. Abell 2029 is unusual in that it hosts 
a strong cooling flow without associated optical line emis- 
sion. This may indicate that the central regions of the cluster 
have been disrupted (perhaps by a minor merger event or 
the strong central radio source in this system), complicat- 
ing the distributions of cooling and absorbing material (cf. 
Section 5.2). 

Four of the CF clusters studied here were also exam- 
ined by White et al. (1991), using Einstein Observatory SSS 
data. These authors measured intrinsic column densities for 
Abell 478, 1795, 2029 and 2142 (in a 3 arcmin radius cir- 
cular aperture, using spectral model C) of 1.7^ h r x 10 21 
atom cm~ 2 , 0.8±g;| x 10 21 atom cm~ 2 , x 10 21 atom 

cm" 2 and l.Sto.l x 10 21 atom cm~ 2 , respectively. (The SSS 
result for Abell 478 has been corrected to account for the dif- 
ferent value of Galactic absorption assumed in that study). 
The intrinsic column densities determined by White et al. 
(1991) are 2 — 5 times larger than the values measured from 
the ASCA data. 

Allen & Fabian (1997) compared the results from their 
X-ray colour deprojection study of PSPC data with the 
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Figure 5. The data, best-fit models and residuals from the fits to the SISO spectra for Abell 1068 and 1795 in the energy range 3.0 — 10.0 
keV, using spectral model A. The best-fit models have then been extrapolated to cover the full 0.6 — 10.0 keV band. The upper panels 
show the data and best-fit models. The lower panels show the residuals to the fits. The figures show that the residuals, which can be 
naturally explained by the introduction of cool emission components with associated intrinsic absorption, take the form of an excess of 
counts between energies of ~ 0.8 — 2.5 keV. 



White el al. (1991) SSS analysis and concluded that if the 
intrinsic absorption were due to cold gas, then agreement 
between the measured column densities could only be ob- 
tained if the covering fraction of the X-ray absorbing mate- 
rial were, in general, < 0.5. We have therefore examined the 
constraints on the covering fraction, /, that can be obtained 
from the ASCA data using spectral model C. The results 
are also listed in Table 9. In all cases we determine a best-fit 
covering fraction for the intrinsic X-ray absorbing material 
of unity. For many clusters, covering fractions of < 70 per 
cent can be firmly ruled out. (The physical significance of 
these results are further explored in Section 6.4). The results 
on the covering fractions, and the reasonable agreement of 
the ASCA and PSPC results for Abell 478 and 1795, suggest 
that the White et al. (1991) results may have systematically 
over-estimated the intrinsic column densities in clusters by 
a factor of a few. 

Before considering in more detail the possible origin and 
nature of the absorbing material, it is pertinent to consider 
whether the intrinsic absorption, inferred to be present using 
a variety of spectral models, could in fact be an artifact due 
to an unconsidered emission process in the clusters. Fig. 
5 shows the SISO spectra for Abell 1068 and 1795, which 
exhibit two of the clearest absorption signatures. The data 
between 3.0 and 10.0 keV have been fitted with a single- 
temperature emission model with Galactic absorption (i.e. 
spectral model A). The residuals to these fits, over the full 
0.6 — 10.0 keV band, are shown in the lower panels. We see 
that the residuals, which can be naturally explained by the 
introduction of cool emission components with associated 
intrinsic absorption (using spectral models C and D), take 
the form of a large excess of counts between ~ 0.8 — 2.5 keV. 
We therefore include the caveat that if some extra emission 



process which can account for such residuals in the X-ray 
spectra were identified, then the requirement for intrinsic 
absorption from the X-ray data would be greatly diminished. 
Finally, we note that the fits with spectral Model A over the 
restricted 3.0 — 10.0 energy range lead to determinations of 
the cluster temperatures for Abell 1068 and 1795 of 6-ll^g 



and 6.1 



+0.5 



keV, respectively, in excellent agreement with 



the results obtained for the full data sets with spectral model 
C. 



6.2 The mass of absorbing material 

The results presented in Section 6.1 show that the cooling 
flows in our sample are typically intrinsically absorbed by 
equivalent hydrogen column densities of a few x 10 21 atom 
cm~ 2 (taking spectral model C as our preferred model). As- 
suming, in the first case, that this absorption is due to cold 
gas with solar metallicity, the mass of absorbing gas, M abs , 
within a radius, r abs , may be approximated as 



M abs ~ 3.2 x 10 7 r2b.AJV H M e 



(2) 



where r a b s is in units of kpc and AAh in 10 2 
For metallicities of OAZq in the absorbing gas (the mean 
emission-weighted value for the X-ray emitting gas in the 
CF clusters; Allen & Fabian 1998b) the implied masses are 
a factor ~ 2 larger. For zero metallicity these masses are a 
factor ~ 3 times larger than the values given by equation 
2. (The scaling factors for different metallicities in the ab- 
sorbing material have determined from spectral simulations 
using the XSPEC code.) 

We have adopted r abs as the radius where the cooling 
time of the cluster gas first exceeds 5 x 10 9 yr in the clusters 
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(a plausible age for the cooling flows in most of the CF sys- 
tems). These radii, and the masses of absorbing gas implied 
by equation 2 (using the column densities determined with 
spectral model C) for both solar metallicity and Z — OAZq, 
are listed in columns 2 — 4 of Table 10. 

The observed masses of absorbing material may be com- 
pared to the masses expected to have been accumulated by 
the cooling flows (M acc ) within the same radii over the past 
5 x 10 9 yr. Assuming that from time t = to i = 5 x 10 9 
yr, the integrated mass deposition rate within radius r a b s 
increases approximately linearly with time, the total mass 
accumulated within radius r a b s after 5 x 10 9 yr is approxi- 
mately 

M acc ~ 2.5M(r < r abs ) x 10 9 M . (3) 

At times t > 5 x 10 9 yr, the mass deposition rate within 
r abs should remain approximately constant. If instead we 
assume M(r < r a b s ) to be constant during the first 5 Gyr, 
the accumulated mass will be twice the value indicated by 
equation 3. (Exactly how the mass deposition rate within 
r a bs grows with time during the first 5 Gyr is unclear and will 
depend upon the evolution of the cluster and the detailed 
properties of the cluster gas.) 

The mass deposition rates within radii r a b s (corrected 
for the effects of absorption due to cold gas, as discussed 
in Section 5.1) and the estimated masses of cooled material 
accumulated by the cooling flows within these radii over a 5 
Gyr period, are summarized in Table 10. Fig. 6 compares the 
accumulated masses with the masses of absorbing material 
determined from the spectral data. The agreement between 
the results is reasonable and is maximized for a metallicity 
in the absorbing gas of ~ 0.4 — O.GZq (depending on the 
growth of the cooling flow over the first 5Gyr) . This supports 
the idea that the observed X-ray absorption may be due to 
material accumulated by the cooling flows. We note that if 
the X-ray absorption in the CF clusters were due to dust 
rather than cold gas, a possibility examined in more detail 
in Section 7, then the agreement between the observed and 
predicted masses of absorbing material shown in Fig. 6 might 
still be expected to hold if the dust were contained in the 
material deposited by the cooling flows. In this case, the 
mass of the absorber calculated with equation 2 would be 
the mass of cooled gas associated with the dust when it was 
deposited from the cooling flow. For a Galactic dust to gas 
ratio, the mass in dust would be ~ 100 times smaller than 
the associated gas mass. 

6.3 The luminosity reprocessed in other 
wavebands 

The luminosities absorbed at X-ray wavelengths must even- 
tually be reprocessed in other wavebands. If the absorbing 
material is dusty (as is likely to be the case in the central 
regions of the clusters e.g. Voit & Donahue 1995; Fabian, 
Johnstone & Daines 1994; Allen et al. 1995) then the bulk 
of this reprocessed emission is likely to emerge in the in- 
frared band {e.g. Dwek, Rephaeli & Mather 1990). Table 11 
summarizes the reprocessed luminosities (i.e. the bolomet- 
ric luminosities absorbed within the clusters) measured with 
spectral model C. The reprocessed luminosities range from 
~ 2 x 10 44 ergs" 1 for Abell 1795 to ~ 5 x 10 45 ergs" 1 for 
RXJ1347.5-1145, with a mean value of ~ 10 45 ergs" 1 . 



Allen et al. (1999) report detections of spatially ex- 
tended 100/mi emission coincident with the X-ray centroids 
in the nearby Centaurus cluster and Abell 2199. The infrared 
fluxes from those clusters were also shown to be in good 
agreement with the predicted values due to reprocessed X- 
ray emission from their cooling flows. We have also used the 
IPAC SCANPI software and archival IRAS data to measure 
the 60 and 100/im fluxes within a four arcmin (radius) aper- 
ture centred on the X-ray centroids for the clusters in the 
present sample. (The medians of the co-added SCANPI re- 
sults were used.) The results are listed in Table 11. The error 
bars associated with the measured fluxes are the root-mean- 
squarc deviations in the residuals, external to the source 
extraction regions, after baseline subtraction. Where no de- 
tection was made, an upper limit equal to three times the 
r.m.s. deviation in the residuals is listed. Where a detection 
was made, we also list the in-scan separations (in arcmin) 
between the peak of the 100/im emission (or the 60/im emis- 
sion in the case of Abell 1835) and the X-ray centre. 

Two of the clusters in our sample, IRAS 09104+4109 
and Abell 1068 have known infrared point sources coincident 
with their X-ray centroids (Kleinmann et al. 1988; Moshir et 
al. 1989). The presence of these sources is confirmed by our 
analysis. The 60 and 100/im data for Abell 1704 and 60/im 
data for Abell 1835 also suggest the presence of unresolved 
sources, coincident with the X-ray centroids for the clusters. 
Abell 478 and Zwicky 3146 have 100/im emission originating 
from close to their X-ray centroids, which appears spatially 
extended, although the infrared flux from Abell 478 is prob- 
ably contaminated by Galactic cirrus. The data for Abell 
2142, 2204, 2261 and 2390 also provide significant detections 
within the four arcmin (radius) source apertures, although 
the peaks of the infrared emission from these systems are 
spatially offset from their X-ray centroids, suggesting that 
the detected flux is likely to originate, at least in part, from 
some other source. 

Several of the clusters included in our study (Abell 586, 
963, 1413, 1795, 2029, 2142 and PKS0745-191) have pre- 
viously been studied at infrared wavelengths, using IRAS 
data, by Wise et al. (1993) and/or Cox, Bregman & 
Schombert (1995). These authors also report no clear de- 
tections of infrared emission from these sources. Edge et al. 
(1999) report measurements of 60pm emission from Abell 
1835 and 60 and 100/im emission from Abell 2390, using a 
similar analysis to that presented here. These authors also 
present detections of 850/mi emission from these clusters 
which, for Abell 1835, they suggest is likely to be due to 
dust heated by the vigorous star formation observed in the 
cluster core or an obscured active galactic nucleus. 

Following Helou et al. (1988) and Wise et al. (1993), we 
can estimate the total infrared luminosities implied by the 
observed IRAS fluxes and flux limits using the relation 

Li-iooo M m ~ 2.8 x 10 44 (-^-) 2 (2.58S 60 + Sioo) ergs" 1 , (4) 

where Seo and Sioo are the 60 and 100/im IRAS fluxes in 
units of Jy. This relation assumes a dust temperature of 
~ 30K (which is consistent with the observations; e.g. Allen 
et al. 1999; Edge et al. 1999) and an emissivity index, n, 
in the range — 2, where the emissivity is proportional to 
the frequency, v n . Where measurements, rather than upper 
limits, were obtained, we associate a systematic uncertainty 
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Figure 6. A comparison of the masses predicted to have been accumulated by the cooling flows within radii r abs over the first 5 Gyr, 
versus the masses within the same radii implied by the observed intrinsic X-ray absorption (assuming the absorption to be due to cold 
gas with a covering fraction of unity). The dashed line indicates equality between the two sets of values. Results are shown both for (left) 
solar metallicity in the absorbing gas and (right) a metallicity of OAZq . 



of 30 per cent with the estimated 1 — 1000/im luminosities, 
which is combined in quadrature with the random errors. 
The 1 — 1000/im luminosities calculated from this relation 
are summarized in Table 11. In general, the measurements 
and upper limits to the observed 1 — lOOO^im luminosities ex- 
ceed the predicted luminosities due to reprocessed X-ray ra- 
diation from the cooling flows. Where measurements (rather 
than upper limits) are made, this may suggest the presence 
of an additional source of infrared flux, such as star forma- 
tion or AGN associated with the central cluster galaxies, 
or some other contaminating (possibly Galactic) source. At 
some level, the massive starbursts at the centres of systems 
like Abell 1068, 1835, 2204 2390 and Zwicky 3146 (Allen 
1995; Crawford et at 1999; Edge et al. 1999) must contribute 
to the detected infrared flux. 



6.4 Systematic uncertainties in the absorption 
measurements 

As discussed in Section 6.1, the measured intrinsic equiv- 
alent hydrogen column densities of absorbing material in- 
ferred from the ASCA spectra exhibit significant system- 
atic variation, depending upon whether the single-phase or 
multiphase emission models are used. We have adopted the 
results obtained with the cooling-flow model (model C) as 
our preferred values, since this model provides a consistent 
description of the spectral and imaging X-ray data for the 
clusters. However, a number of systematic uncertainties af- 
fecting the absorption results remain with the cooling flow 
model. 

Spectral model C undoubtedly over-simplifies the situa- 
tion in a real cooling flow. In particular, the model assumes 
that the absorbing material lies in a uniform screen in front 



of the cooling flows. If the absorber is intrinsically associ- 
ated with the clusters, then the absorbing material must, 
to some extent, be distributed throughout the cooling flows. 
Allen & Fabian (1997; see also Wise & Sarazin 1999) discuss 
two different, limiting geometries for the absorbing material; 
partial covering and multilayer absorption. The partial cov- 
ering model, with a covering factor / < 1, can be expected 
to apply or where there are < 1 individual absorbing clouds 
along each line of sight and/or where the absorbing material 
is strongly clumped on large scales. However, as discussed in 
Section 6.1, covering fractions significantly less than unity 
are found to provide a poor description of the ASCA spec- 
tra, suggesting that the absorbing material is likely to be 
more smoothly distributed throughout the cooling flows. 

The multilayer absorption model discussed by Allen & 
Fabian (1997) is applicable where the absorbing material 
is made up of a large number of similarly sized absorb- 
ing clouds, each with a column density much less than the 
total intrinsic column density (AJVh), homogeneously dis- 
tributed throughout the X-ray emitting medium. For an 
intrinsic emission spectrum A(E), the observed spectrum, 
A'(E), emerging from the emitting/absorbing region may 
be written (where cr(E) is the absorption cross-section) as 



A'(E) = A(E)( 



e -a(E)AN H \ 

V a(E)AN H ) ' 



(5) 



We have simulated cooling flow spectra with multilayer 
absorption, with total column densities of between 10 21 and 
5 x 10 22 atom cm~ 2 , and fitted these with simple models 
which assume that the absorbing material acts as a uniform 
screen in front of the cooling flow (as with model C) . For the 
simulations we adopt an emission spectrum appropriate for a 
constant pressure cooling flow, with an upper temperature of 
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7 keV and a metallicity of OAZq. We use a response matrix 
appropriate for the ASCA SIS detectors. The comparison 
between the fitted (uniform screen) column densities and 
the true total (multilayer) values are shown in Fig. 7. We see 
that for multilayer column densities of 10 21 — 5 x 10 22 atom 
cm~ 2 , the fitted column densities underestimate the true 
values by factors of 5 — 10. (These results are not sensitive 
to any reasonable choice of Galactic column density.) 

If the multilayer model approximates the true distri- 
bution of X-ray emitting and absorbing material in cooling 
flows, then the results obtained with spectral model C would 
imply true equivalent column-densities of intrinsic absorb- 
ing gas in the cluster cores of, typically, a few 10 22 atom 
cm~ . Importantly, we also note that for signal-to-noise ra- 
tios appropriate for ASCA observations of bright clusters 
(with count rates of ~ let s _1 in the SIS detectors and expo- 
sure times of ~ 40 ks) , the simple uniform-screen absorption 
model provides a good fit to the simulated spectra with mul- 
tilayer absorption, for total column densities < 3 x 10 22 atom 
cm~ 2 . (For larger total column densities, complex residuals 
are detected in the spectra at energies < 2.0keV.) The large 
intrinsic column densities implied by the multilayer models 
would have strong implications for the nature of the absorb- 
ing matter. In particular, column densities of > 10 22 atom 
cm~ 2 of gaseous material with solar metallicity would be 
very difficult to reconcile with current HI and CO limits, 
even if this material were very cold and highly molecular. 
(We also note that in the case of a multilayer distribution, 
the implied masses of absorbing gas would often exceed the 
predicted accumulated masses due to the cooling flows and 
imply, at least in part, some other origin for the absorbing 
medium.) In the following section, we explore these issues 
in more detail and summarize the current observational con- 
straints on the physical nature of the absorbing matter. 



7 SEARCHES FOR THE X-RAY ABSORBING 
MATERIAL IN OTHER WAVEBANDS 

The results presented in this paper provide a consistent pic- 
ture for the X-ray properties of cooling flows in X-ray lumi- 
nous clusters of galaxies. The ASCA spectra and ROSAT im- 
ages provide independent determinations of the mass cool- 
ing rates in the clusters in good agreement with each other, 
when a consistent method of analysis is employed (Section 
5.1). The spectral data require the presence of large col- 
umn densities of intrinsic X-ray absorbing material associ- 
ated with the cooling gas (Section 6). The implied masses 
of absorbing material are in reasonable agreement with the 
masses expected to have been accumulated by the cooling 
flows over their lifetimes (Section 6.2), identifying this ma- 
terial as a plausible sink for the bulk of the cooled gas from 
the flows. Despite the evidence for intrinsic X-ray absorption 
associated with cooling flows, however, searches for this ma- 
terial in other wavebands have, to date, proved largely un- 
successful (at least outside of the central 5 — 20 kpc). We 
here briefly review these observations and comment on their 
implications for the physical state of the X-ray absorbing 
matter. More detailed discussions of a number of these issues 
are given by Ferland, Fabian & Johnstone (1994), Daines, 
Fabian & Thomas (1994), Fabian et al. (1994a), O'Dea et 
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Figure 7. The apparent column density determined from fits to 
cooling flow spectra with a uniform-screen absorption model (e.g. 
spectral model C) when in reality the absorber is distributed in 
a multilayer manner (Allen & Fabian 1997). The vertical axis 
shows the fitted column density and the horizontal axis the true 
total (multilayer) column density through the emitting/absorbing 
region. Both axes are in units of 10 21 atom cm -2 . 



al. (1994a), Voit & Donahue (1995), O'Dea & Baum (1996), 
Arnaud & Mushotzky (1988) and Henkel & Wiklind (1998). 

As mentioned in Section 5.2, the central (5 — 20 kpc) 
regions of CF clusters are commonly observed to be sites 
of ongoing star formation with powerful, associated opti- 
cal line emission (e.g. Johnstone et al. 1987; Heckman et 
al. 1989; McNamara & O'Connell 1989; Crawford & Fabian 
1992; Allen 1995; Cardiel et al. 1995, 1998; McNamara et al. 
1996; Voit & Donahue 1997; Crawford et al. 1999). Nearly all 
CF clusters with short central cooling times (t coo i i; 2Gyr) 
exhibit such phenomena (Peres et al. 1998), whereas sys- 
tems with longer central cooling times (including all NCF 
clusters) generally do not. The observed levels of star for- 
mation imply the presence of at least moderate masses of 
cooled gas in the cores of CF clusters, sufficient in a few 
cases to account for the bulk of the mass deposited by the 
cooling flows in those regions (Allen 1995) and at least some 
causal connection with the cooling flows. 

At radio wavelengths, extensive searches have been car- 
ried out for 21cm line-emission associated with warm atomic 
hydrogen in the cores of CF clusters (e.g. Valentijn & Gio- 
vanelli 1982; O'Dea & Baum 1996). The negative results 
obtained (with the exception of M87; Jaffe 1992) typically 
constrain the column densities of optically thin material to 
be less than a few 10 19 atom cm -2 and, for an assumed spin 
temperature, Ts ~ 20K, the number of clouds along the line 
of sight to be < 1 (O'Dea & Baum 1996). 

Searches have also been made for the 21 cm absorp- 
tion signature of cold atomic hydrogen (e.g. McNamara, 
Bregman & O'Connell 1990; Jaffe 1992; Dwarakanath, van 
Gorkom & Owen 1994; O'Dea, Gallimore & Baum 1995; 
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O'Dea & Baum 1996; Johnstone et al. 1998). Such mate- 
rial has been detected across the central lew kpc ol a lew 
CF clusters, including Abell 426 (Crane, Van der Hulst & 
Haschick 1982; Sijbring et al. 1989; Jaffe 1990) 2A0335+096 
(McNamara et al. 1990), Hydra A (Dwarakanath, van 
Gorkom & Owen 1995), MKW3s (McNamara et al. 1990) 
and Abell 2597 (O'Dea, Baum & Gallimore 1994) with typi- 
cal column densities ol a lew 10 20 (Ts/100) atom cm~ 2 . How- 
ever, lor most clusters only negative results have been ob- 
tained, with upper limits to the column densities ol 10K gas 
(assuming a covering fraction ol unity and a velocity width 
of ~ 1.5 kms -1 ; Dwarakanath et al. 1994, O'Dea & Baum 
1996) ol 10 18 - 10 19 atom cm" 2 . 

Laor (1997) presents tight limits on the column den- 
sity of atomic hydrogen towards the core of Abell 426 
(AJVh <4x 10 17 atom cm~ 2 ) from Lya observations made 
with the Hubble Space Telescope (see also Johnstone & 
Fabian 1995). Together, the results on HI emission and ab- 
sorption, and Lya absorption, imply that the bulk ol the 
material responsible for the observed X-ray absorption can- 
not be in the form of atomic hydrogen. 

An issue of crucial importance in interpreting the 21cm 
HI results and the results from CO studies (see below) is the 
expected temperature ol X-ray absorber, il in gaseous form, 
which has proved controversial. Johnstone, Fabian & Tay- 
lor (1998; updating earlier calculations by Ferland, Fabian 
& Johnstone 1994 and Fabian et al. 1994a) suggest that 
cooled gas clouds, in equilibrium with the cooling-flow envi- 
ronment, are likely to be very cold. Johnstone et al. (1998) 
suggest that dust-free clouds form an extended outer enve- 
lope, with a column density ol ~ 5 x 10 21 atom cm~ 2 and a 
temperature of 13-17K, beyond which the gas drops to the 
microwave background temperature. For Galactic dust/gas 
ratios, however, the results are significantly modified, such 
that the cloud temperature drops to the microwave back- 
ground temperature after a column density of only ~ 10 20 
atom cm~ . At such low temperatures, dusty gas is likely 
to be highly molecular. However, these predictions differ 
from those of O'Dea et al. (1994a) and Voit & Donahue 
(1995) who suggest a minimum temperature for such clouds 
of ~ 20K, even in the presence of dust. Braine et al. (1995) 
also suggest a minimum temperature of ~ 10K for the ab- 
sorbing gas. Possible reasons for the discrepancies between 
these results (with contrasting views) are discussed by Voit 
& Donahue (1995) and Ferland et al. (in preparation). 

Extensive searches have been carried out for CO as- 
sociated with molecular gas in cooling flows. To date CO 
emission has only been detected in Abell 426 (Lazareff et 
al. 1989; Mirabel, Sanders D.B. & Kazes I. 1989). Braine et 
al. (1995) also present limits on CO absorption in the cen- 
tral regions of this cluster). The upper limits to the column 
density of molecular hydrogen in the inner regions of other 
clusters (for an assumed kinetic temperature of Tq Q — 20K 
and a covering fraction of unity) are typically < a few 10 20 
atom cm -2 (O'Dea et al. 1994a, Antonucci & Barvainis 1994, 
McNamara & Jaffe 1994; Braine & Dupraz 1995; O'Dea & 
Baum 1996) which is significantly below the X-ray inferred 
column densities. The CO results also constrain the number 
of clouds along the line of sight (for t£ = 20K) to be < 10 
(O'Dea et al. 1994a; O'Dea & Baum 1996). The CO emission 
limits appear to exclude the possibility of large column den- 
sities {i.e. values consistent with the X-ray measurements) 



of molecular gas, with a kinetic temperature > 10K. How- 
ever, if the bulk of the X-ray gas cools to the microwave 
background temperature, as suggested by Johnstone et al. 
(1998), the CO limits may be consistent with the X-ray data. 

It should be noted that the determinations ol molecu- 
lar hydrogen masses from CO observations are affected by a 
number of additional uncertainties. Firstly, the CO/H2 ra- 
tios in cluster cooling flows may differ from Galactic values. 
In particular, where the temperature and metallicity of the 
absorbing gas is low, the use of a standard Galactic CO/H 2 
ratio may underestimate the mass ol molecular gas (Maloney 
& Black 1988; Madden et al. 1997). 

It is unclear whether X-ray absorbing gas accumulated 
by cooling flows could have substantially sub-solar metal- 
licity. Allen & Fabian (1998b) determine a mean emission- 
weighted metallicity for the X-ray emitting gas in CF clus- 
ters of Z ~ 0.4Zq . These authors also show that CF clusters 
generally contain metallicity gradients, so that the metal- 
licity of the X-ray gas within the cooling radii of CF sys- 
tems will, on average, exceed this value. Thus, it seems un- 
likely that the material currently being deposited from the 
cooling flows will have substantially sub-solar metallicities. 
However, if the X-ray absorbing material is formed from the 
matter deposited by the cooling flows over their lifetimes, 
the metallicity ol this material could be somewhat lower. 
Within an inhomogeneous cooling flow {e.g. Nulsen 1986, 
1998; Thomas, Fabian & Nulsen 1987) the denser gas at 
any particular radius will tend to have a smaller volume 
filling factor. If the metals are evenly distributed through- 
out the volume of the cooling flow at any particular radius 
(whether this will occur is unclear) , the least dense gas will 
contain most of the metals (Reisenegger, Miralda-Escude & 
Waxman 1996). Since the densest material in cooling flows 
will have the shortest cooling time and be deposited first, 
the material accumulated by the cooling flows over their his- 
tories could have a significantly lower metallicity than the 
material being deposited today. Conversely, if the metals 
in cooling flows were concentrated in the densest material, 
the metallicity of the X-ray absorbing material could exceed 
that of the X-ray emitting gas. 

Fabian et al. (1994a), have suggested that the masses 
ol molecular gas inferred from CO observations could be 
significantly underestimated if most of the CO in the ab- 
sorbing material has frozen onto dust grains. Such freezing 
should occur rapidly (on timescales ~ 10 5 yr) if even small 
amounts of dust are initially present in the absorbing gas 
(Voit & Donahue 1995). Daines et al. (1994) and Fabian et 
al. (1994a) discuss how dust is likely to form even in clouds 
that are initially free of dust, on timescales ~ 10 9 yr (al- 
though see also Voit & Donahue 1995). Optical, UV and 
sub-mm studies ol the cores ol cooling flows show that dust 
lanes {e.g. Sparks, Macchetto & Golombek 1989; McNamara 
& O'Connell 1992) and intrinsic reddening (Hu 1992; Allen 
et al. 1995; Edge et al. 1999) are common. 

Elston & Maloney (1994) and Jaffe & Bremer (1997) 
report detections of H2(1-0)S(1) emission from warm (~ 
2000K) molecular hydrogen from K-band spectroscopy of 
the inner few kpc of a number of CF systems (including 
Abell 478 and and PKS0745-191). No detections of such 
emission have been made in NCF systems. The H2S emission 
is likely to arise from the heated skin of cold molecular clouds 
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Figure 8. (a) The equivalent column densities of hot hydrogen ions through two simulated clusters, one CF and one NCF system. The 
gas distributions have been modeled as isothermal spheres, with core radii of 50 kpc (CF; upper curve) and 300 kpc (NCF; lower curve) 
respectively, and have been normalized to provide total X-ray gas masses within radii of 3Mpc of 2 X 1O 14 M0. The distribution of 
X-ray emitting gas is much more centrally concentrated in the CF cluster, particularly within the central 100 — 200 kpc. (b) The X-ray 
absorbing column density, as a function of radius, due to dust distributed throughout the simulated CF (upper curve) and NCF (lower 
curve) clusters. The dust is assumed to have a multilayer distribution within the X-ray emitting gas, but be modeled as a simple uniform 
screen in front of the clusters. 



which could also be responsible for the optical emission-line 
phenomena discussed above. 

A final possibility, first discussed in detail by Voit & 
Donahue (1995; see also Arnaud & Mushotzky 1998), is that 
the X-ray absorption could be due to dust, with little or 
no associated gas. The K-absorption edge of oxygen is a 
primary source of X-ray absorption and its signature is fairly 
insensitive to the physical state of the oxygen, so long as 
it is not highly ionized. The introduction of a simple OIK 
absorption edge at E ~ 0.54 keV in the ASCA analysis 
(such as might be associated with oxygen-rich, silicate dust 
grains) typically provides at least as good a model for the 
intrinsic absorption as a cold gaseous absorber. Arnaud & 
Mushotzky (1998) have shown that 0.35 - 7.0 keV Broad 
Band X-ray Telescope data for the Perseus cluster require 
excess absorption which is significantly better-explained by 
a simple oxygen absorption edge (which they associate with 
dust) than by a cold, gaseous absorption model. Verification 
of this important result will be possible in the near future 
using observations made with the Chandra Observatory and 
XMM. 

Optical and UV spectroscopy of the central emission- 
line nebulosities in cooling flows indicate the presence of 
significant amounts of dust (Hu 1992; Allen et al. 1995; 
Crawford et al. 1999). Allen et al. (1995) also show that 
the column densities of absorbing gas inferred from the op- 
tical reddening studies are in reasonable agreement with the 
X-ray values, for Galactic dust/gas ratios. The infrared ob- 
servations reported here and elsewhere, together with the 
sub-mm results of Edge et al. (1999), require the presence 
of significant dust masses in the core regions of at least a 



few nearby and/or exceptionally massive cooling flows (for 
which the best data exist) and are consistent with dust be- 
ing a common feature of CF clusters. Limits on background 
counts of galaxies and quasars behind clusters are consistent 
with extinctions ranging from Ab = 0.2 — 0.5 (e.g. Boyle, 
Fong & Shanks 1988; Romani & Maoz 1992), although mea- 
surements of galaxy colours in nearby clusters constrain the 
reddening on cluster- wide scales to be E(B — V) < 0.06 mag 
(Ferguson 1993). 

Voit & Donahue (1995) argue that the levels of dust 
required to account for the observed X-ray absorption are 
unlikely to lie entirely within the cores of the clusters, and 
should be distributed more widely throughout the cluster 
gas. ROSAT results on the spatial distribution of the X- 
ray absorbing material (Allen & Fabian 1997) show it to be 
centrally concentrated within the cooling radii of clusters 
and possibly be confined within cooling flows. However, the 
X-ray data are also consistent with a scenario in which the 
absorbing material is distributed in a manner similar to the 
X-ray emitting gas. This is illustrated in Fig. 8(a), where 
we show the equivalent column densities of hydrogen ions in 
the X-ray emitting gas through two simulated clusters; one 
CF and one NCF system. The gas distributions have been 
modeled as isothermal spheres, with core radii of 50 and 
300 kpc, respectively (and have been normalized to provide 
total X-ray gas masses within 3Mpc of the cluster centres of 
2 x 10 14 Mq). We see that the distribution of X-ray emitting 
gas is much more centrally concentrated in the CF system, 
particularly within the central ~ 100 kpc. 

Fig. 8(b) shows the corresponding X-ray absorbing col- 
umn densities as a function of radius (as would be deter- 
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mined from observations with ASCA-like CCD detectors in 
the 0.6 — 10.0 keV band) in the case where the intraclus- 
ter gas has an equivalent metallicity in dust grains of ~ 0.2 
solar. For simplicity, we assume that the dust has an absorp- 
tion spectrum in the ASCA band similar to that of metal- 
rich cold gas (Morrison & McCammon 1983; qualitatively 
similar results are obtained for an oxygen edge at 0.54 keV). 
We also assume that the absorber follows a multilayer dis- 
tribution within the X-ray emitting medium, but is mod- 
eled as a uniform absorbing screen in front of the emitting 
regions. The results imply mean emission- weighted column 
densities (where we assume the X-ray emissivity to be pro- 
portional to square of the gas density) of ~ 2 x 10 21 atom 
cm~ 2 for the CF cluster, and a few xlO 20 atom cm~ 2 for 
the NCF system. (The absolute results on the excess column 
densities will depend upon the chemical composition of the 
dust). We thus see that if the X-ray absorption is due to 
dust, with a spatial distribution that follows the X-ray gas, 
then we can expect absorption signatures to be present in 
both CF and NCF clusters, but to be stronger in CF systems 
(assuming in each case that an appropriate spectral model 
is used in the analysis). We would also expect the absorp- 
tion to be concentrated towards the cores of CF clusters, in 
agreement with the observations. We note that if the X-ray 
absorption were entirely due to a pre-existing distribution 
of dust grains distributed throughout the clusters, then the 
agreement between the observed masses of absorbing matter 
(calculated assuming that the absorption is due to cold gas 
confined within the cooling flows) and the masses predicted 
to have been accumulated by the cooling flows over their 
lifetimes (Fig. 6) must be regarded as coincidence. However, 
it remains possible that both material accumulated by the 
cooling flows and a more extended distribution of large dust 
grains (presumably due to supernovae enrichment at early 
epochs) could contribute to the X-ray absorption. 

The lifetime of grains of radius a/im to sputtering in hot 
gas of density n is ~ 2 x 10 a/n yr (Draine & Salpeter 1979). 
Provided that individual grains exceed 10/im in radius, they 
should survive for a Hubble time or longer throughout NCF 
clusters and beyond the cooling radius in clusters with cool- 
ing flows. Within cooling flows, the gas density rises inward 
so the grains are increasingly sputtered, releasing the metals 
into the gas phase which thus becomes increasingly metal 
rich towards the cluster centre. This provides one possi- 
ble explanation for the abundance gradients inferred to be 
present in most CF clusters (Allen & Fabian 1998b; see also 
Irwin & Bregman 1999). However, 10/im grains are optically 
thick to soft X-rays and so some distribution of grain sizes 
extending to smaller values is required if dust is responsi- 
ble for the observed X-ray absorption {e.g. Laor & Draine 
1993). We note that the sputtering of small (a < 0.1/im) 
grains will significantly modify the reddening law in clusters 
with respect to the standard Galactic relation, reducing the 
optical/UV reddening (Laor & Draine 1993). 

In conclusion, the available data from other wavebands 
provide some support for the large column densities of intrin- 
sic absorbing material inferred to be present from the X-ray 
data, at least in the innermost regions of cooling flows. Opti- 
cal and sub- mm observations of star formation in the cores of 
CF clusters and K-band observations of emission lines from 
molecular hydrogen suggest that significant masses of dusty, 
molecular gas are present in the central few tens of kpc of 



many CF clusters. The constraints from 21cm observations 
suggest that the bulk of the material accumulated by cool- 
ing flows cannot remain in a long-lived reservoir of atomic 
hydrogen. Substantial masses of molecular gas may be dis- 
tributed throughout cooling flows and as yet have avoided 
detection, but this material must be very cold (T ~ 3K) 
and dusty, a possibility that remains controversial. The very 
large X-ray column densities (> 10 22 atom cm~ 2 ) required 
if the absorbing material has a multilayer distribution are 
probably incompatible with a gaseous absorber. It remains 
possible that dust grains, either present in the cooled mate- 
rial deposited by the cooling flows or distributed throughout 
the clusters in a manner similar to the X-ray gas, may be 
responsible for much of the observed X-ray absorption. 

8 CONCLUSIONS 

The main conclusions that may be drawn from this paper 
may be summarized as follows: 

(i) We have demonstrated the need for multiphase mod- 
els to consistently explain the spectral and imaging X-ray 
data for the CF clusters included in our study. The mass de- 
position rates from the cooling flows, independently inferred 
from multiphase analyses of the ASCA spectra and depro- 
jection analyses of the ROSAT HRI images, exhibit good 
agreement, especially once the effects of intrinsic X-ray ab- 
sorption have been accounted for in a consistent manner. 
The mass deposition rates from the largest cooling flows ex- 
ceed 1000 M yr -1 , identifying these as some of the most 
massive cooling flows known. 

(ii) We have confirmed the presence of intrinsic X-ray 
absorption in the cluster spectra using a variety of spectral 
models. The measured equivalent hydrogen column densi- 
ties are sensitive to the spectral models used in the analysis, 
ranging from a few 10 20 atom cm~ 2 for a simple isothermal 
emission model to a few 10 21 atom cm -2 using our preferred 
cooling-flow models, assuming in each case that the absorber 
lies in a uniform foreground screen. Both the CF and NCF 
systems exhibit excess absorption at a similar level (on aver- 
age) when analysed with the same, simple isothermal model. 

(iii) The masses of X-ray absorbing material inferred 
to be present in the CF clusters (assuming the absorption 
to be due to cold gas with a covering fraction of unity) are 
in reasonable agreement with the masses expected to have 
been accumulated by the cooling flows over their lifetimes. 

(iv) The ASCA spectra constrain the covering fraction 
of the absorbing material acting on the cooling flows to be 
close to (or exceed) unity. If the X-ray absorption is due to 
many small, similarly-sized clouds along each line of sight, 
intermixed with the X-ray emitting gas, then the column 
densities inferred from the spectral analysis may signifi- 
cantly underestimate the true column densities of absorbing 
material in the cooling flows. 

(vii) We have summarized the constraints on the phys- 
ical properties of the X-ray absorbing material from ob- 
servations in other wavebands. Substantial 60 and 100/im 
fluxes, sufficient to account for the bolometric luminosities 
absorbed and reprocessed within the clusters, are detected 
from several of the largest cooling flows. Optical observa- 
tions of star formation in the cores of cooling flows also 
suggests that significant masses of molecular gas are present 
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in the central few tens of kpc of many CF systems. Con- 
straints from 21cm observations suggest that the bulk of 
the material accumulated within the cooling radii cannot re- 
main in a long-lived reservoir of atomic hydrogen. Substan- 
tial masses of molecular gas may be distributed throughout 
cooling flows and as yet have avoided detection, although 
this material must be very cold (T ~ 3K), a possibility that 
remains controversial. The very large X-ray column densities 
(> io 22 atom cm -2 ) required if the absorbing material has 
a multilayer distribution are probably incompatible with a 
gaseous absorber, ft remains possible that dust grains, either 
present in the material accumulated by the cooling flows or 
distributed more widely in a manner similar to the X-ray 
gas, may be responsible for much of the X-ray absorption. 
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Table 1. Summary of the ASCA Observations 



Cluster 


z 


iV H 


Date 


SO 


SI 


G2 


G3 


Abell 2744 


0.308 


0.16 


1994 Jul 04 


37605 


26086 


62749 


62753 


Abell 478+ 


0.088 


3.00 


1994 Feb 24 


27953 




34560 


34554 


Abell 520 


0.203 


0.78 


1994 Sep 09 


16802 


15871 


17656 


17617 


Abell 586 


0.171 


0.52 


1994 Mar 22 


8482 




17668 


17648 


PKS0745-191 


0.103 


4.24 


1993 Nov 06 


29146 




37553 


37553 


Abell 665 


0.182 


0.42 


1993 Sep 18 


22926 


19869 


34815 


34799 


IRAS 09104+4109 


0.442 


0.10 


1993 Nov 12 


33447 


24348 


38655 


38655 


Abell 773 


0.217 


0.14 


1994 Nov 22 


38500 


38086 


40841 


40831 


Abell 963 


0.206 


0.14 


1993 Apr 22 


29611 


29039 


29883 


29881 


Zwicky 3146 


0.291 


0.30 


1993 May 18 


29663 


28643 


32622 


32632 


Abell 1068 


0.139 


0.14 


1996 Dec 06 


23777 


23645 


21463 


21639 


Abell 1413 


0.143 


0.22 


1993 Dec 11 


26084 


19862 


36173 


36175 


Abell 1689 


0.184 


0.18 


1993 Jun 26 


29575 


23642 


37817 


37817 


Abell 1704 


0.216 


0.18 


1994 Apr 16 


18847 


11784 


20598 


20596 


RXJ1347.5-1145 


0.451 


0.49 


1995 Jan 17 


27882 


17549 


38968 


38958 


Abell 1795 


0.063 


0.12 


1993 Jun 16 


31284 




37649 


37641 


MS1358. 4+6245 


0.327 


0.19 


1995 Apr 27 


32532 


30815 


31981 


31513 


Abell 1835 #1 


0.252 


0.23 


1994 Jul 20 


18051 


17532 


17464 


17460 


Abell 1835 #2 






1994 Jul 21 


16876 


16444 


16412 


16410 


MS1455.0+2232 


0.258 


0.32 


1994 Jul 18 


29388 


28539 


28903 


28895 


Abell 2029 


0.077 


0.31 


1994 Feb 19 


34211 


32420 


35038 


34963 


Abell 2142 


0.089 


0.42 


1994 Feb 21 


12335 




15768 


15766 


Abell 2163 


0.208 


1.21 


1993 Aug 08 


25126 


18224 


32760 


32322 


Abell 2204 


0.152 


0.57 


1994 Aug 20 


12749 


11865 


14698 


14698 


Abell 2218 


0.175 


0.32 


1993 Apr 30 


28241 


26054 


37970 


37968 


Abell 2219 


0.228 


0.18 


1994 Aug 07 


32705 


31697 


35849 


35849 


Abell 2261 


0.224 


0.33 


1996 Feb 24 


19224 


18823 


19249 


19219 


Abell 2319 


0.056 


0.80 


1993 Jul 21 


9842 


5007 


13169 


13404 


MS2137.3-2353 


0.313 


0.36 


1994 May 08 


15167 


15732 


17035 


17056 


Abell 2390 


0.233 


0.68 


1994 Nov 13 


6172 


2632 


10340 


10338 


AC114 


0.312 


0.13 


1995 Nov 09 


36739 


36295 


35987 


35971 



Notes: Columns 2 and 3 list the redshift and Galactic column density, in units of 10 21 atom cm -2 from Dickey & Lockman (1990), with 
the exception of Abell 478+, for which the value is from Allen & Fabian (1997). For PKS0745-191, this study (Section 3.2) suggests a 
Galactic column density of ~ 3.5 X 10 21 atom cm -2 , which is lower than the Dickey & Lockman (1990) value. Column 4 lists the dates 
of the ASCA observation. Columns 5 — 8 list the net exposure times (in seconds) in each of the four ASCA detectors, after all screening 
and cleaning procedures were carried out. 
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Table 2. Extraction radii for the ASCA spectra 



Cluster SO SI Chip mode 

(amin/kpc) (amin/kpc) 



AKr.ll 97 44 


Q 4/11 AC\ 

0.4:/ 1 14:U 


9 7 /Qft^ 
z. ( / yuo 


4C1 "1 
1(1) 


AKoll A7S 


O.o/ QUO 






AKoll ^90 
/\DC11 OZU 


Q /Q99 


Q 1 /7Q/L 
O. 1/ / y4: 


9^'9^ 

V / 


.MDtJll ooo 


Q /CI A 

O.O/ 014: 




4C1 "1 




4:. z/ uoo 








4^/1 n7n 

4:.0/ 1U 1 U 


O.O/ 1ZOU 




TRA^ (1Q1 044-41 OQ 


4 0/1 fi90 

4:.U/ 1UZU 


< -±..\j 1 1UZU 


4(4) 


AKpII 773 
.MUCH I f o 


4 4 / 1 1 sn 

L ±. L ± 1 lOU 


o.y/ iuou 




AHpII 


3.6/932 


9 8 /79 c i 

Z. O/ 1 ZU 






O.O/ 1 1 zu 


4 3 /I 3Q0 
4.0/ ioyu 


2(2) 


Ahf>ll lOfik 

iT_Ut_-ll 1UUO 


0.4:/ UOO 


9 fi /^09 
z. u/ ouz 




AKoll 1413 

rt_[JCll 14:10 


4 fl /7Qfl 
4.0/ i yu 


4 fl /7Q0 
4.0/ i yu 


4(4) 


Aholl IfiSQ 


4 O/Q 1 ^ 
i. o/ you 


i. u/ yoo 


4(4) 


Abell 1704 


3.4/909 


2.6/695 


4( , 11 


RXJ1347.5-1145 


4.0/1630 


4.0/1630 


4(1) 


Abcll 1795 


4.5/445 




4(2) 


MS1358.4+6245 


3.9/1350 


2.9/1010 


4(2) 


Abell 1835 #1 


4.1/1220 


3.2/948 


2(1) 


Abell 1835 #2 


4.1/1220 


3.2/948 


2(1) 


MS1455.0+2232 


3.7/1110 


2.8/842 


2(1) 


Abcll 2029 


3.5/413 


2.7/319 


2(1) 


Abell 2142 


6.0/802 




4(4) 


Abell 2163 


3.7/964 


4.5/1170 


4(2) 


Abell 2204 


4.9/1010 


4.7/973 


2(1) 


Abell 2218 


3.7/852 


2.5/575 


4(1) 


Abell 2219 


5.3/1470 


4.5/1250 


2(1) 


Abell 2261 


4.5/1230 


3.6/987 


2(2) 


Abell 2319 


4.2/373 


4.0/355 


4(1) 


MS2137.3-2353 


4.3/1460 


3.4/1150 


2(2) 


Abell 2390 


3.6/1010 


2.5/704 


4(1) 


AC114 


3.4/1150 


2.9/980 


1 



The radii of the circular extraction regions for the ASCA SIS data 
(in arcmin and kpc) and the chip modes used in the observations 
(cither 1,2 or 4-CCD mode). The numbers in parentheses indicate 
the number of chips contributing to the extracted spectra. For the 
CIS data a fixed extraction radius of 6 arcmin was used. For Abcll 
2142, the 2 arcmin (3 arcmin) radius region surrounding the X-ray 
bright Seyfert-1 galaxy 1556+274 was masked out and excluded from 
the analysis for the SIS (CIS) data. 
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Table 3. Summary of the ROSAT Observations 



Cluster 


Date 


Exposure (s) 


RA. (J2000). 


Dec. (J2000. 


Abell 2744 


1994 Dec 09 


34256 


00 h 14 m 18.7 s 


-30°23'11" 


Abell 478 


1991 Feb 10 


22840 


04 h 13 m 25.4 s 


10° 27' 58" 


Abell 520 


1994 Mar 05 


12728 


04 h 54 m 10.1 s 


02°55'27" 


Abell 586 


1994 Oct 06 


23136 


07 h 32 m 20.4 s 


31°37'55" 


PKS0745-191 


1992 Oct 20 


23750 


07 h 47 m 31.1 s 


-19° 17' 47" 


Abell 665(P) 


1991 Apr 10 


38308 


08 h 30 m 58.9 s 


65°50'37" 


IRAS09104+4109 #1 


1994 Nov 08 


7968 


09 h 13 m 45.5 s 


40° 56' 29" 


IRAS09104+4109 #2 


1995 Apr 12 


21904 






Abell 773 


1995 Apr 15 


6240 


09 h 17 m 53.4 s 


51°43'29" 


Abell 963 


1992 Nov 24 


10104 


10 h 17 m 03.4 s 


39°02'51" 


Zwicky 3146 #1 


1992 Nov 27 


15214 


10 h 23 m 39.8 s 


04°11'11" 


Zwicky 3146 #2 


1993 May 17 


10831 






Abell 1068 


1995 Apr 29 


15760 


10 h 40 m 44.6 s 


39° 57' 12" 


Abell 1413(P) 


1991 Nov 27 


7696 


ll h 55 m 18.7 s 


23° 24' 12" 


Abell 1689 #1 


1994 Jul 22 


13080 


13 h ll m 29.1 s 


-01°20'40" 


Abell 1689 #2 


1995 Jun 24 


9648 






Abell 1704 


1995 Apr 17 


41360 


13 h 14 m 24.8 a 


64°34'39" 


RXJ1347.5-1145 


1995 Jan 28 


15760 


13 h 47 m 31 s 


-11°45'11" 


Abell 1795 #1 


1992 Jun 25 


2768 


13 h 48 m 52.7 s 


26°35'27" 


Abell 1795 #2 


1993 Jan 21 


11088 






Abell 1795 #3 


1994 Jun 23 


11080 






MS1358. 4+6245 #1 


1991 Nov 05 


2528 


13 h 59 m 51.0 s 


62°31'04" 


MS1358. 4+6245 #2 


1993 May 14 


15872 






Abell 1835 


1993 Jan 22 


2850 


14 h 01 m 02.0 s 


02°52'40" 


MS1455. 0+2232 #1 


1992 Jan 11 


4088 


14 h 57 m 15.0 s 


22°20'36" 


MS1455. 0+2232 #2 


1993 Jan 20 


4230 






MS1455. 0+2232 #3 


1994 Jul 07 


6584 






Abell 2029 


1990 Jul 25 


17758 


15 h 10 m 56.2 s 


05°44'42" 


Abell 2142 


1994 Jul 26 


19776 


15 h 58 m 20.1 s 


27°13'52" 


A 1 11 o ~1 /"> O 

Abell 2163 


1994 Aug 13 


36248 


16 15 45.9 


— Oo 08 58 


Abell 2204 


1995 Sep 01 


15472 


16 h 32 m 47.0 s 


05°34'33" 


Abell 2218 #1 


1994 Jan 05 


11520 


16 h 35 m 52.5 s 


66° 12' 29" 


Abell 2218 #2 


1994 Apr 05 


27016 


77 77 


77 77 


Abell 2218 #3 


1994 Jun 14 


30016 


77 77 


77 77 


Abell 2218 #4 


1994 Jun 17 


24304 


77 77 


77 77 


Abell 2219 


1994 Jan 17 


13242 


16 h 40 m 20.2 s 


46°42'29" 


Abell 2261 


1995 Feb 11 


16288 


17 h 22 m 27.3 s 


32°07'58" 


Abell 2319 


1991 Mar 22 


5444 


19 h 21 m 12.2 s 


43°56'38" 


MS2137.3-2353 


1994 Apr 24 


13656 


21 h 40 m 15.2 s 


-23°39'41" 


Abell 2390 


1993 Nov 23 


27764 


21 h 53 m 36.5 s 


17°41'45" 


AC114 #1 


1993 May 17 


10048 


22 h 58 m 48.7 s 


-34°48'19" 


AC114 #2 


1994 May 09 


13144 


77 77 


77 77 



Notes: Columns 2 and 3 list the dates of observation and exposure times (in seconds). Where more than a single observation of a cluster 
was made, details for each observation are listed. Columns 4 and 5 give the coordinates of the centroids of the X-ray emission from the 
clusters. For Abell 665 and 1413, HRI images were not available at the time of writing and PSPC data were used for the imaging analysis. 
This is indicated by a (P) after the cluster name. 
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Table 4. Results from the spectral analysis 



Parameters 


Model A 


Model B 


Model C 


Model D 


kT 
Z 

Abell 2744 N H 


11 04+ - 79 
ir.ut_ 73 

n 1V+0.09 
u,± ' -0.09 

0.16 


7 ?c .+0.59 
'• '°-0.53 

16+ 006 

1 44+0.21 




o c+0.9 

».t>_ 7 

n 17 +0.08 
u - 1( -0.07 
i oc+0.14 
± - zo -0.22 


M S 










kT 2 








< 1.7 


AiV H 








> 5.0 


L x 


31.3 


31.0 




30.9 


LboI 


70.1 


67.8 




U.C. 


X 2 /DOF 


862.7/737 


739.2/736 




729.0/733 



kT 6.02±°-li 6.09l°i« 8.1±£"* 7.13±°|| 

Z 0.32+n 2 3 0.321°^ °- 35 io':o4 0-34 + °;^ 

Abell 478 AT H 3.00 2.94+° Jj* 3.00 3.00 

M s 13471^ 

kT 2 ' 0.76+°' ^ 

AJV H 1.7±°"jj 9.4+Jj 

L x 24.2 24.2 24.2 24.2 

L B oi 53.8 53.5 63.6 71.6 

X 2 /DOF 914.3/815 912.8/814 879.5/813 842.0/812 



Abell 520 



kT 
Z 
N H 


o Qo+0.76 
°- oo -0.67 

14+0 10 

"■^-O.IO 

0.78 


7 94+0-96 
1 -0.90 

n QO+0.23 
u - a -0.21 






M s 










kT 2 










AN H 










Lx 


15.8 


15.7 








34.1 


33.9 






X 2 /DOF 


431.3/430 


430.2/429 







Abell 586 



kT 


7 irc+0.72 
'• ±o -0.62 


7 02+°- 94 

1 - u -0.80 


10.7+ 3 ° 9 3 


> 7.1 


Z 


31+0 13 


Sl+°' 13 
U-01_o .12 


o.36±g-;i« 


S7+°- 17 
u - cw -0.16 


N H 


0.52 


n KO+0 35 
u - oy -0.32 


0.52 


0.52 


M S 






U.C. 




kT 2 








9+ 4 ' 3 
u.a_ 2 


AN n 






U.C. 


1 q+67.8 


Lx 


10.7 


10.7 


10.7 


10.7 


LboI 


22.7 


22.7 


26.1 


31.8 


X 2 /DOF 


207.6/232 


207.4/231 


205.2/230 


202.4/229 



PKS0745-191 



kT 


6 7l+°- 14 

D -' i -0.14 


6 47+°- 17 


8.7+1;° 


7 O+0.7 


Z 


n ni+0.03 


O1+0.03 
u " 31 -0.02 


U -" JO _0.03 


S6+ 04 
U " 3O -0.04 


N H 


3.50 


o fjn+0.10 


o q 7 +0.46 
d-*' -0.51 


3 fil +0 22 
• 3 - Di -0.18 


M S 










kT 2 








i n+0.5 
L " i -0.3 


AN H 








8-0 + 3 5 5 ° 


Lx 


29.5 


29.4 


29.5 


29.6 


LboI 


65.2 


65.3 


74.0 


77.6 


X 2 /DOF 


1326/1235 


1315/1234 


1296/1232 


1270/1231 
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Table 4. Spectral Results - continued 



Parameters 



Model A 



Model B 



Model C 



Model D 





kT 


q no + 0.58 
y - UJ _0.52 


o 19 +0.62 
8-12_ 54 


— 


14.2+ 5 ° 2 4 




Z 


22 +ao7 


n 90+0. 06 
u -^°_0 07 


— 


0.331°^ 


Abcll 665 


N H 


0.42 


72+° :i4 


— 


o.o±°* 




M S 


— 





— 


— 




kT 2 


— 


— 


— 


< 4.4 




AN K 


— 


— 


— 


> 2.5 




Lx 


17.8 


17.6 




17.7 




LboI 


38.5 


37.9 


— 


U.C. 




v 2 /DOF 


753.1/745 


740.4/744 




732.4/741 




kT 


6-2S±0 ; 5| 


in an 

°- OM -0.62 


S -°-1.8 


9.3+J 7 




Z 




44+ ' 12 
u -^-0.12 


n c-i +0.17 


0-67l o j 2 


IRAS 09104+4109 


N H 


0.10 


S1+ ' 26 


0.10 


0.10 




M s 


— 




1655l 6 7 tl 


— 




kT 2 














AN H 


— 


— 


3-0±l;I 






Lx 


22.0 


21.9 


21.8 


21.9 




LboI 


46.3 


46.5 


55.1 


65.3 




X /DOl 


401.5/361 


399.5/360 


392.7/359 


382.6/358 




kT 


q 9q +0.69 
y -^ y _0 60 




— 


— 




Z 


0- 21+ 009 
— u.uy 


0- 21+ 008 

— U.Uo 


— 


— 


Abell 773 


N H 


0.14 


0.42+ /" 
— 0.14 








Ma 












kT 2 












AN H 


— 


— 


— 


— 




Lx 


14.8 


14.5 










32.0 


31.1 


— 


— 




x /DUb 


736.0/648 


724.8/647 








kT 






6-13l° ; 6 6o" 


°- 1J _0.30 




Z 


0.301°;°? 


0.30ig-°? 


0.31±g;S 


so+° 11 


Abell 963 


N H 


0.14 


16 +0 14 


0.14 


0.14 




M s 


— 





U.C. 


— 




kT 2 


— 





— 


U.C. 




AN H 


— 


— 


U.C. 


U.C. 




Lx 


12.7 


12.7 


12.7 


12.7 




LboI 


27.0 


27.0 


27.0 


27.0 




X 2 /DOF 


566.6/544 


566.5/543 


566.6/542 


566.6/541 




kT 


6 4l+°- 26 


6-8o±8:ii 


11.3+2 7 


10.6+1° 




Z 


n 97 +0.06 
u - z < -0.06 


n 9 7 +0.06 


Q on + 0.07 


n o 7 +0.08 
u -' 3 ' -0.09 


Zwicky 3146 


N H 


0.30 


H+ - 11 
u - ±J --o.io 


0.30 


0.30 




M s 






22281^7 






kT 2 








2 6+ 1 ' 




AN H 






-1 9+0.3 


4+ 1 ' 6 




Lx 


36.9 


37.2 


37.3 


37.4 




LboI 


78.6 


78.8 


91.8 


88.7 




X 2 /DOF 


724.1/735 


716.3/734 


699.9/733 


696.0/732 
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Table 4. Spectral Results - continued 





Parameters 


Model A 


Model B 


Model C 


Model D 


Abell 1068 


kT 
Z 
N H 
M S 
kT 2 
AN H 
Lx 


4 - zi -0.14 
43+0-08 

0.14 
7.00 


3 46+ 16 
49+O.O8 

99+° 15 

u - y -0.14 

6.66 


r = + 1.4 
°-°-0.9 

42+ - 10 
-0.08 

0.14 

1345li« 

, O+l.O 

•^-o.s 
6.80 


16.6+8*6 

0.14 

9 C-+0.3 
Z -°-0.6 

1 7+ 1 - 1 

1 -'-0.4 

6.95 




LboI 


16.1 


17.1 


21.8 


20.4 




X 2 /DOF 


728.4/552 


613.6/551 


607.6/550 


585.1/549 



Abell 1413 



kT 


7 54+0.29 
' -° -0.27 


6 64+ - 30 

°-° -0.28 


8 5+ 13 
°-°-0.8 


9.8±tl 


Z 


n 90+O.O5 
U - z8 -0.05 


n 90+O.O5 
U - z8 -0.05 


S0+ 06 
U - 3U -0.05 


0.36+°;^ 


N H 


0.22 


66+° 10 

u -° -0.11 


0.22 


0.22 


M S 






644+^ 




kT 2 








94+l.S 
Z ' -0.8 


AN H 






34 +i.3 

°-^-0.8 


q q+3.6 


Lx 


15.5 


15.3 


15.4 


15.4 




32.9 


32.8 


36.4 


36.6 


X 2 /DOF 


992.4/954 


938.9/953 


934.3/952 


924.7/951 



Abell 1689 


kT 
Z 
N H 
M S 
kT 2 
AN H 


9 66+ ' 37 

29+0' 05 
u - zy _0.05 

0.18 


1O+0-44 
y - ly -0.39 
n 90+O.O5 
u ' zy -0.05 
n 09+O.O9 
^•"-0.08 


10.01+ 2 
30+ 06 

0.18 
,r-i+291 
001 -209 

4 1+5-6 


10 4+ ' 8 
34+ 06 
0.18 




Lx 


32.4 


32.3 


32.2 


32.2 




LboI 


70.3 


69.6 


62.6 


84.5 




X 2 /DOF 


1085/1125 


1077/1124 


1075/1123 


1065/1122 



Abell 1704 



kT 


4 7 ,+0.38 


4 20+ - 45 
^• zu _0.40 


r- 7 +3.5 
°-'-1.3 


6 1+ 6 ' 8 
D.±_ 1 2 


Z 


34+0- 15 


30+O.16 
u ' oy -0.13 


n ,0+0.17 
0.38_ 16 


39+ 14 
u - oa -0.17 


N H 


0.18 


64+ ' 31 

u -° -0.29 


0.18 


0.18 


M S 






9S7+ 382 
y °'-616 




kT 2 








1 0+ 1 ' 8 


AN H 






, ,+5.2 
•'•°-1.6 


10.9+^ 


Lx 


8.32 


8.13 


8.19 


8.28 


LboI 


18.6 


19.0 


22.9 


30.5 


X 2 /DOF 


249.1/281 


241.9/280 


240.2/279 


234.5/278 



RXJ1347.5-1145 



kT 


12.46lg;|g 


10.44+ ,' 88 


26.4+^ 2 8 3 


19 q+24.5 


Z 


0.38±g;iS 


0-34±°;g| 


Q 43+0.11 


55+ 20 
u - oo -0.17 


N H 


0.49 


0.98+^ 


0.49 


0.49 


M S 






3479+335 
°^' 9 -1151 




kT 2 








2.7+ 2 ; 2 


AN H 






9 7+3.O 
Z -'-0.7 


8.3+ 5 3 3 2 


Lx 


94.4 


93.9 


93.5 


93.5 


LboI 


217.8 


210.7 


225.0 


261.8 


X 2 /DOF 


803.0/741 


772.7/740 


764.8/739 


763.9/738 
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Table 4. Spectral Results - continued 





Parameters 


Model A 


Model B 


Model C 


Model D 


Abell 1795 


kT 
Z 

N H 
M S 
kT 2 
AN H 
Lx 


5 40+ 08 

n sfi+°' 03 

U-O0_ .02 
0.12 

11.0 


C qq+0.10 

36+ 03 
n 17+0-05 

U ' 1 ' -0.05 

11.0 


r- o 7 +0.29 
°-°'-0.25 
36+ 03 

0.12 

300t»> 4 

n i +0.6 
11.0 


7.i± :. 

0.381^4 
0.12 

2 5+ 1 ' 3 
7 +1 '° 
10.40 




LboI 


23.6 


23.6 


25.2 


25.6 




X 2 /DOF 


1442/1236 


1440/1235 


1422/1234 


1388/1233 



MS1358.4+6245 



kT 


7 48+0-83 


r 7^+0.74 




U.C. 


Z 


S2+° 15 


35+° ' 14 


0.38±°;1« 


0.36t°£ 


N H 


0.19 


1 04+ ' 32 
1-U -0.29 


0.19 


0.19 


M S 






69ll 3 f 7 




kT 2 








5 8+ ' 7 
°'°-0.6 


AN H 






6 4+8-7 

D -*-3.8 


1 8+ 1 ' 9 


Lx 


10.8 


10.5 


10.5 


10.6 


LboI 


22.8 


22.9 


26.3 


24.4 


X 2 /DOF 


333.1/321 


308.1/320 


311.3/319 


305.1/318 



Abell 1835 



kT 


o 91+O.3I 


7 qq + 0.35 

' -"-o.so 


9-8±i:l 


11 1+ 4 - 4 


Z 


oc+0.06 


34+ 05 


40+ 06 


47+0 08 

u -^< -0.08 


Nh 


0.23 


u - oo -0.09 


0.23 


0.23 


M S 






1761± 5 23 




kT 2 








2 9+ 1 ' 2 
z ' a -0.9 


AN H 






q 9+1.6 
^■ Z -0.8 


o 1+3.1 
• i - L -1.6 


Lx 


45.4 


44.8 


44.9 


44.9 


LboI 


96.3 


95.1 


106.2 


105.8 


X 2 /DOF 


1292/1311 


1255/1310 


1247/1309 


1239/1308 



MS1455.0+2232 



kT 


4 83+ ' 22 


4 qo+0.27 


5 6+ 3 - 1 


7 n +25.0 
'' u -2.0 


Z 


n 09+0. 09 
U - JZ -0.07 


n qk+0.08 
°" 36 -0.07 


n o 7 +0.09 
u -°'-0.08 


30+0.10 


Nh 


0.32 


73+0 17 


0.32 


0.32 


M S 






1728t™ 




kT 2 










AN H 








1.81- 


Lx 


16.6 


16.3 


16.4 


16.7 


LboI 


37.0 


37.8 


44.7 


39.8 


X 2 /DOF 


489.7/495 


472.7/494 


471.3/493 


468.8/492 



Abell 2029 



kT 


7 7 9+0.13 
'•' z -0.13 




8 47+0- 41 
-0.36 


8.95+° j? 


Z 


43+0 03 


0.42lH 3 3 


46+ 03 
u -™-0.03 




Nh 


0.31 


58+ 04 
U - 08 -0.05 


0.31 


0.31 


M S 






547« 2 




kT 2 








i-6±g:i 


AN H 






q 1 +0.5 


6.6l 2 ;l 


Lx 


19.9 


19.9 


19.9 


19.9 


LboI 


41.8 


42.0 


45.5 


47.1 


X 2 /DOF 


1401/1068 


1289/1067 


1256/1066 


1226/1065 
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Table 4. Spectral Results - continued 



Parameters 



Model A 



Model B 



Model C 



Model D 





kT 


a fi7 +0.32 


<-? OO+0.35 
'•° 2 -0.32 


a.J_ 7 






Z 


25+ ' 05 
u - z0 — 0.05 




n 9 7 +0.05 
u - z ' —0.05 


0.30±g;g? 


Abell 2142 




0.42 


°- 77+ o :i io 

— U. 1U 


0.42 


0.42 




M s 






529+lIT 
J J — ±QA 






kT 2 


— 


— 




9+ 1 ' 6 

u -0.3 




AN H 






3 6+ 1 ' 8 


ir n + 7.5 
i0 - U -11.4 




Lx 


30.4 


30.4 


30.4 


30.4 




LboI 


64.6 


64.6 


67.8 


80.1 




X 2 /DOF 


737.5/649 


701.2/648 


696.5/647 


684.1/646 




kT 


13.83+°-™ 


10.85l°;S 




11.811-° 




Z 


23+ ' 07 


n 99 +0.06 
u - zz _0.06 




0.23l°;» 7 


Abell 2163 


N H 


1.21 


i 07+O.I3 




1 76+°' 13 

1 - '"-0.14 




M S 












kT 2 








< 1.5 




AN H 








> 10 




L x 


61.0 


60.1 




60.2 




LboI 


147.0 


137.3 




U.C. 




X 2 /DOF 


1138/1068 


1047/1067 




1037/1064 





kT 


7 4n+ - 30 

'• 4U -0.28 


6 23+ 30 

"■^-0.28 


9.2+H 


13.7 +21 7 3 




Z 


41+ 06 

u -0.06 


40+ 06 

u '^ u -0.05 


46 +o:or 

u -™-0.07 


"■ o,i -0.09 


Abell 2204 


N H 


0.57 


1 1S+ 012 


0.57 


0.57 




M s 






2103l 3 7 « 






kT 2 








3 - 8 -a9 




AN H 






o O+0.9 

°- J -0.8 


1 9+ 1 ' 4 




Lx 


35.2 


34.5 


34.6 


34. 7 




LboI 


74.2 


73.4 


85.5 


82.9 




X 2 /DOF 


1004/834 


921.7/833 


916.3/832 


909.8/831 




kT 


7 n5 +0.36 
'• UO -0.35 


7 10+0. 50 
'■18-0.45 








Z 




0-18+°;^ 






Abell 2218 


N H 


0.32 


o.26l° ;il 








M S 












kT 2 












AiV H 












Lx 


10.8 


10.8 








LboI 


23.1 


23.1 








X 2 /DOF 


710.0/694 


709.4/693 








kT 


12 42+ - 77 
-^■^ -0.69 


9 46+ ' 63 
y - 4D -0.57 




17 9 +29.2 




Z 


n 10+0. 08 


20+ ' 05 
u - zu _0.06 




26+° 15 
U ' ZD -0.11 


Abell 2219 


N H 


0.18 


77+O.IO 
u - ' ' -0.10 




41+ 031 
u -^ 1 -0.41 




M S 












kT 2 








1 6+ 1 ' 4 

1 -°-0.6 




AN H 








7 4+4-7 
'■ 4 -4.2 




Lx 


38.9 


38.0 




38.3 




LboI 


89.4 


83.3 




110.8 




X 2 /DOF 


991.0/898 


892.6/897 




861.9/894 
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Table 4. Spectral Results - continued 



Parameters 



Model A 



Model B 



Model C 



Model D 





kT 


8.82±°;£ 


6 64+ ' 51 


10.9t 5 2 9 2 


13.4+1 D 




Z 




S1+ 07 
u - Ji -0.08 


n o 7 +0.09 


52+ ' 17 


Abell 2261 


N H 


0.33 


1 17+0-17 
11 ' -0.17 


0.33 


0.33 




M S 


— 




1323^263 


— 




kT 2 


— 


— 


— 


Z '°-0.7 




AJV H 






3.8^-j 






Lx 


23.9 


23.2 


23.3 


23.3 






51.3 


50.0 


58.3 


61.2 




X 2 /DOF 


604.8/550 


522.0/549 


515.5/548 


512.2/547 




kT 


9 30+ ' 41 
a - ou -0.38 


8 57+ - 45 




n-3±i;I 




Z 


n OO+0.06 


n 00+O.O6 
u - oo -0.06 




n o 7 +0.09 
U - J '-0.08 


Abell 2319 


N H 


0.80 


1 05+° 11 




f) QQ+ - 23 




M S 












kT 2 








noi+0.47 
u - ai -0.27 




AJV H 












L x 


16.9 


16.9 




16.9 




LboI 


36.8 


36.7 




43.4 




X 2 /DOF 


1156/1086 


1140/1085 




1106/1082 



MS2137.3-2353 



kT 


5 1 6+0-39 
°- lo -0.34 


4.4ilS:S 


5.2lo:? 


U.C. 


Z 


44+°' 14 


0.48±g;i| 


0.50±g;i5 


55+0- 24 


N H 


0.36 


n qt;+0.29 


0.36 


0.36 


M s 






1467t 8 7 ^ 




kT 2 










AJV H 








O q+47.7 

• 3 - J -2.5 


Lx 


17.0 


16.6 


16.6 


16.7 


LboI 


36.9 


38.1 


44.1 


40.9 


X 2 /DOF 


239.5/279 


226.1/278 


226.3/277 


223.5/276 





kT 


10.13±J; 22 


8 OK+l-30 
°- ao -1.04 


14.5+l 5 2 5 


16.0+f 4 ° 




Z 


o sfi+o- 15 

u - ao -0.16 


34+0- 14 


o.4oi° ;« 


55+0- 34 
U.OO_ Q 24 


Abell 2390 


N H 


0.68 


-i 11+0.32 
- L -- L1 -0.30 


0.68 


0.68 




M S 






1 e;'M+ 581 

iooo_ 1112 






kT 2 








9 c+3.2 
Z -°-1.4 




AN H 






2 9+ 7 ' 6 


5 5+ 14 ' 9 
°-°-4.0 




Lx 


41.3 


41.0 


41.0 


41.1 




LboI 


90.8 


89.2 


101.1 


109.2 




X 2 /DOF 


327.5/344 


321.9/343 


320.4/342 


318.3/341 




kT 


9 76+ 1 ' 04 

a -' u -0.85 










Z 


20+ ' 12 
u - zu -0.13 


20+ 10 






AC114 


N H 


0.13 


64+ ' 22 

-0.21 








M S 












kT 2 












AiV H 












Lx 


17.5 


17.2 








LboI 


38.1 


37.0 








X 2 /DOF 


491.1/436 


474.4/435 
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Table 4. Spectral Results 



The best-fit parameter values and 90 per cent (Ax 2 = 2.71) confidence limits from the spectral analysis of the ASCA data. Temperatures 
(fcT), metallicities (Z), column densities (AT H ) and intrinsic column densities (ANn) were linked to take the same values in all four 
detectors. With spectral models C and D, the mass deposition rates of the cooling flows (Ms) an d luminosities of the cooler spectral 
components were also linked to take the same values, but scaled by a normalization factor proportional to the total flux measured in 
the detector (Section 3.1). Only the normalization of the hotter isothermal emission component was allowed to vary independently for 
each detector. Temperatures are in keV, metallicities as a fraction of the solar value (Anders & Grevesse 1989), column densities and 
intrinsic column densities in units of 10 21 atom cm~ 2 and mass deposition rates in Mq yr — 1 . The Ms values for model C are measured 
with the G3 detector. Lx values are the 2 — 10 keV X-ray luminosities in the rest frame of the source. LboI values are the bolomctric 
luminosities, corrected for the effects of Galactic and intrinsic absorption. All luminosities are in units of 10 44 ergs _1 and are measured 
in the G3 detector. An entry of U.C. indicates that a parameter was unconstrained by the data. Note that the analysis of PKS0745-191 
with spectral model C has the Galactic column density included as a free parameter (Section 3.2). The fit to Abell 2390 with spectral 
Model D has a second possible solution with kT ~ 0.2 keV and kT2 ~ 9.0 keV which provides a comparable x 2 > although this solution 
appears physically implausible. 
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Table 5. The goodness of fit for the different spectral models 
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GOODNESS OF FIT IMPROVEMENT 



Cluster 


Model A 


Model B 


Model C 


Model D 


A— > B 


B^ CF 


B— > 2 r 


Abell 2744 


9.9 x 10~ 4 


0.46 


— 


— 


1.00 


— 


0.98 


Abell 478 


8.6 x 10~ 3 


8.8 x 10" 3 


5.2 x 10~ 2 


0.23 


0.77 


1.00 


1.00 


Abell 520 


0.47 


0.48 


— 


— 


0.69 


— 


0.57 


Abell 586 


0.87 


0.87 


0.88 


0.90 


0.34 


0.74 


0.90 


PKS0745-191 


1.4 x 10~ 3 


5.4 x 10~ 2 


0.10 


0.22 


1.00 


1.00 


1.00 


Abell 665 


0.43 


0.53 


— 


— 


1.00 


— 


0.96 


IRAS 09104+4109 


7.0 x 10~ 2 


7.4 x 10~ 2 


0.11 


0.18 


0.83 


0.96 


1.00 


Abell 773 


9.1 x 10~ 3 


1.8 x 10~ 2 


— 


— 


1.00 


— 


0.82 


Abell 963 


0.24 


0.24 


0.23 


0.22 


0.24 


0.00 


0.00 


Zwicky 3146 


0.61 


0.67 


0.81 


0.83 


0.99 


1.00 


1.00 


Abell 1068 


6.0 x 10~ 7 


3.3 x 10~ 2 


4.5 x 10~ 2 


0.14 


1.00 


0.93 


1.00 


Abell 1413 


0.19 


0.62 


0.65 


0.72 


1.00 


0.91 


1.00 


Abell 1689 


0.80 


0.84 


0.85 


0.89 


1.00 


0.84 


0.99 


Abell 1704 


0.92 


0.95 


0.96 


0.97 


0.90 


0.63 


0.97 


RXJ1347.5-1145 


5.7 x 10~ 2 


0.20 


0.25 


0.26 


1.00 


0.98 


0.96 


Abell 1795 


4.0 x 10~ 5 


4.3 x 10" B 


1.5 x 10~ 4 


1.3 x 10~ 3 


0.83 


1.00 


1.00 


MS1358.4+6245 


0.31 


0.67 


0.61 


0.69 


1.00 


0.00 


0.64 


Abell 1835 


0.64 


0.86 


0.89 


0.91 


1.00 


0.99 


1.00 


MS1455.0+2232 


0.56 


0.75 


0.75 


0.77 


1.00 


0.52 


0.80 


Abell 2029 


< io- 7 


3.0 x 10" 6 


4.7 x 10~ 5 


4.2 x 10~ 4 


1.00 


1.00 


1.00 


Abell 2142 


8.9 x 10~ 3 


7.3 x 10~ 2 


8.7 x 10~ 2 


0.15 


1.00 


0.92 


1.00 


Abell 2163 


6.7 x 10~ 2 


0.66 






1.00 




0.98 


Abell 2204 


4.3 x 10" 5 


1.7 x 10~ 2 


2.2 x 10~ 2 


2.9 x 10~ 2 


1.00 


0.97 


0.99 


Abell 2218 


0.33 


0.33 






0.55 




0.50 


Abell 2219 


1.6 x 10~ 2 


0.54 






1.00 




1.00 


Abell 2261 


5.3 x IO" 2 


0.79 


0.84 


0.85 


1.00 


0.97 


0.99 


Abell 2319 


6.9 x 10~ 2 


0.12 






1.00 




1.00 


MS2137.3-2353 


0.96 


0.99 


0.99 


0.99 


1.00 


0.00 


0.65 


Abell 2390 


0.73 


0.79 


0.79 


0.81 


0.97 


0.55 


0.88 


AC114 


3.5 x IO" 2 


9.3 x 10~ 2 






1.00 




0.20 



The goodness of fit results for spectral models A— D, listing the probabilities of exceeding the best-fitting reduced \ 2 values, assuming 
in each case that the model correctly describes the data. Column 6 lists the significance of the improvements to the fits obtained with 
model B over model A {i.e. including the Galactic column density as a free parameter in the isothermal model). Columns 7 and 8 list 
the significance of the improvements obtained by using the multiphase spectral models [cither the cooling flow (CF) or two-temperature 
models (2T)] in preference to the single-phase model B. 
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Table 6. Results from the deprojection analysis of the ROSAT images 



MASS PROFILES DEPROJECTION RESULTS 

Cluster binsizc a r c *cool *100 r cool Mi 

(arcscc/kpc) (kms -1 ) (kpc) (10 9 yr) (10 9 yr) (kpc) (Mgyr 1 ) 



Abell 2744 


16/89.4 


Abell 478 


12/26.5 


Abell 520 


24/102 


Abell 586 


16/60.3 


PKS0745-191 


12/30.3 


Abell 665 


30/118 


IRAS 09104+4109 


12/81.0 


Abell 773 


16/71.5 


Abell 963 


16/69.0 


Zwicky 3146 


12/64.8 


Abell 1068 


12/38.6 


Abell 1413 


30/98.7 


Abell 1689 


12/47.7 


Abell 1704 


12/53.5 


RXJ1347.5-1145 


12/81.7 


Abell 1795 


16/26.3 


MS1358.4+6245 


8/46.3 


Abell 1835 


8/39.5 


MS1455.0+2232 


8/40.1 


Abell 2029 


16/31.5 


Abell 2142 


16/35.7 


Abell 2163 


12/52.1 


Abell 2204 


12/41.4 


Abell 2218 


16/61.4 


Abell 2219 


10/46.2 


Abell 2261 


8/36.6 


Abell 2319 


32/47.4 


MS2137.3-2353 


12/67.7 


Abell 2390 


16/75.1 


AC 114 


12/67.6 



11201*° 


450 


18.8lf 4 ° 


8601™ 


70 




9801™ 


450 


16.7l 2 6 +° 


1050l«° 


130 


5 5+°' 7 
°-°-0.7 


930l 9 ° 


37.5 


1 l+°+ 


930l 3 3 ° 


200 




1120±So 


45 


2.01°;! 


9901™ 


180 


9.8l« 6 7 


750lf 


80 


4 1+1-3 


1150t»2 


50 


i 7 +o.i 
1 -'-o.i 


780+?° 


30 


l-21 :l 


9601™ 


130 


8.6!°;? 


9901^ 


80 


2 g+0.5 


750lf ° 


30 


2 1+°+ 


ISSOio'™ 


75 


2 6+°' 2 
Z ' D -0.1 


7401^ 


50 


1 4+0-1 

!'^-0.1 


830l 3 8u ° 


40 


9 «+l- 2 


1000 + J,o 


50 


i £+0.5 
l-°-0.3 


760+ 7 ?,° 


45 


1 l+°- 2 
I- 1 -!).! 


880+*° 


80 


1 5+° + 


9101™ 


100 


4 3+°' 8 


1210l 3 4 ° 


300 


19 7 +12.3 
lz -'-4.6 


1020l 6 3 ° 


25 


o.94l°;°l 


830l 3 ° 


230 


10.41*+ 


10501™ 


250 


9-OlS 8 


1030l 2 « 


80 


O Q+1.4 


890+2°, 


150 


a c + 10.4 

y.o_ 3 2 


8301^° 


32 


1 2+°+ 
l- z -().l 


1190lg° 


60 


4 2+°' 3 
^• z -0.2 


10301°° 


300 


16 Q+ 32 - 7 
lo.»_ 6 2 



19.2 


< 134 


< 101 


3.11 


999+30 
^ -50 


81011™ 


16.3 


< 154 


< 85 


7.94 
2.65 


119t£ 


159+°° 
650li 2 ° 


13.0 


< 177 


< 232 


2.32 


1931?° 


620llf 


11.2 


< 179 


< 238 


5.60 
2.53 
2.75 
8.63 
4.34 
3.59 
3.00 
3.42 
4.76 
2.44 
2.15 
3.54 
5.38 


188+i 3 

— 84 

205l 2 3 
185+^ 

— 50 
110+29 
119 _70 

180lgj 
187+? 3 

— 54 

196+?° 
1881JS 

1 ,0 + 187 
llo_ 44 

231+ 2 ° 

193l 2 2 
154!^ 


340+l 2 ° 

— 183 

10571- 7 
408l 9 28 
137lf 5 
5631^5 

306l? 2 
1378l™ 6 

449lf 6 

10 o+369 
1AS -50 
H54+432 
1101_482 

712+g 5 
576+?° 
303l 7 f 


12.1 


< 130 


< 90 


2.01 


235lg 


1007+ 98 
iuu ' -263 


12.8 


< 153 


< 133 


10.7 


< 254 


< 485 


5.63 
9.81 
1.65 
5.45 


1251™ 

1° 9 ^39 
216+ 4 f 

1461^ 


227+ g? 8 
801" 
754+ 2 °? 
3701°^ 


16.6 


< 135 


< 120 



Notes: Column 2 lists the binsize (in arcsec and kpc) used in the analysis. Columns 3 and 4 summarize the velocity dispersions and core 
radii used to parameterize the mass profiles (Section 4). Columns 5 and 6 list the cooling times i.e. the time for the gas to cool from 
the ambient cluster temperature at constant pressure (in units of 10 9 yr) within the central bin and within a fixed radius of 100 kpc 
(*cool an d <i00i respectively). Cooling radii (r coo i) are the radii where the cooling time of the cluster gas becomes equal to the age of 
the Universe (1.3 X 10 10 yr). Mass deposition rates (.Mi) are the integrated mass deposition rates within the cooling radii in Mq yr~ x . 
Errors on the velocity dispersions are 90 per cent confidence limits. Errors on the cooling times are the 10 and 90 percentile values from 
100 Monte Carlo simulations, using the best-fit mass models for the clusters. The upper and lower confidence limits on the cooling radii 
are the points where the 10 and 90 percentiles exceed and become less than the Hubble time, respectively. Errors on the mass deposition 
rates are the 90 and 10 percentile values at the upper and lower limits for the cooling radius. For Abell 478, 586, 1795, 2029 and 2219, 
the introduction of a second mass component was found to improve the isothcrmality of the temperature profile and for these clusters the 
multi-component mass models (described in Section 4) have been used in determining the deprojection results. We note that the errors 
on the deprojection results do not account for the uncertainties in the total mass profiles. 
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Table 7. Absorption-corrected mass deposition rates 





INTRINSIC 


ORIGINAL 


CORRECTED 


SPECTRAL 




ABSORPTION 


DEPROJ. 


DEPROJ. 


ANALYSIS 


Cluster 


A7V H 


Mi 


M c 


M S 




(10 21 atom cm" 2 ) 


(Moyr- 1 ) 


(Meyr 1 ) 


(M yr- 1 ) 



Abell 478 






1059t 42 J 


1347l 2 i 7 


PKS0745-191 






9501200 


14551^6 


IRAS 09104+4109 


-LI fi 

°- u -1.3 


4-1 zL7 

620±Jf 


1060+222 


1655l 7 « 
2228lgr 


Zwicky 3146 


i O+0.3 
1 -^-0.3 


1057+ 117 


1358t 449 5 


Abell 1068 


q o+l.l) 
a "*-0.8 






1345+i 4 7 


Abell 1413 


3 4+1.3 


137+^ 


32lti? 


644+^° 


Abell 1689 


4 1+5-6 
4t - x — 1.8 




20161^5 


nr., +291 


Abell 1704 


q Q + 5.2 

• 3 -°-1.6 


306t 9 7 2 


6321™ 


987t 3 6 f 6 


RXJ1347.5-1145 


9 7+3.0 
^•'-0.7 


1378±™ 6 


17901 4 ; 3 


3479+f 3 ^ 


Abell 1795 


9 1+0.6 


449lf 6 


10381 229 


300+ 1 ? 4 


MS1358.4+6245 


6 4+8-7 

D - 4 -3.8 


19 o+369 
iZJ -50 


7n q+1668 
' u ' 3 _445 


691+ 34 ? 


Abell 1835 


q 9+1-6 


1154+ 432 


21H+2154 
zl±i -921 


1761±g 3 


MS1455.0+2232 


4 1+ 43 

4t - x — 1.6 


712+ 53 5 


1 227+ 715 
-172 


17281S 4 


Abell 2029 


q 1+0.5 


576±» 


1220+. 223 


547+ 82 

— 81 


Abell 2142 


3 fi+!-8 


303+. 7 f 


1303t 4 °l 


529+ 177 


Abell 2204 


q q+0.9 


1007+^3 


1660lg 4 


2103l 35 7 6 8 


Abell 2261 




227+ 238 


680+ 2 f 4 


1323l 2 « 


MS2137.3-2353 


6-2^i 


754+263 
'0^_117 


1537lJ 4 f 


1467± 880 , 


Abell 2390 


2.9+TJ 




5qq+ 852 
Jaa -163 


1533l»l 2 



A comparison of the mass deposition rates measured from the deprojection analysis, corrected for the effects of intrinsic absorption due 
to cold gas (Mc), with the values determined from the spectral analysis (Ms)- Errors on the corrected deprojection values account for 
the statistical uncertainties in both the deprojected quantities and intrinsic column densities. No entries arc listed for Abell 586 or 963 
since the spectral data do not constrain the mass deposition rates and intrinsic column densities for these systems. 
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Table 8. Cooling flow flux fractions 



Cluster (L coo i/Lx) 



+0.15 
-0.12 
+0.19 
-0.17 
+0.45 
-0.32 
+0.37 
-0.31 
+0.26 
-0.20 
+0.13 
-0.09 
+0.10 
-0.06 
+0.52 
-0.36 
+0.10 
-0.57 
+0.05 
-0.05 
+0.62 
-0.18 
+0.22 
-0.15 
+0.60 
-0.24 
+0.04 
-0.04 
+0.07 
-0.05 
+0.25 
-0.14 
+0.18 
-0.20 
+0.40 
-0.16 
+0.49 
-0.42 



AKpII 47K 


0.37 


PK"^074^ 1 Q1 


0.36 




0.55' 




0.63 


Abcll 1f)fiS 


0.74" 


A Koll 1 /1 1 Q 
/VDCll 1410 


n qiy 
u.ou 


Abell 1689 


o.io; 


Abcll 1704 


0.48; 


RXJ1347.5-1145 


o.9o; 


Abcll 1795 


o.n; 


MS1358.4+6245 


0.38; 


Abcll 1835 


0.34; 


MS1455.0+2232 


0.40; 


Abell 2029 


0.21; 


Abell 2142 


0.13; 


Abell 2204 


0.49; 


Abcll 2261 


0.59; 


MS2137.3-2353 


0.29; 


Abcll 2390 


0.51; 



Notes: The fractions of the 2 — 10 kcV X-ray luminosities contributed 
by the cooling flows (determined with spectral model C and measured 
in the G3 detector). The errors are 90 per cent confidence limits. 
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Table 9. The column densities of intrinsic X-ray absorbing material 



CF CLUSTERS NCFs 
Cluster B-A C-A C'-A D-A / Cluster B-A 



Abell 478 


n nc+0.08 







QO+0.28 

° -0.38 


9 4+1-6 
a - 4 -1.4 


> 0.85 


Abell 2744 


i 90+0.21 

1 - zS -0.21 


Abell 586 


07+ ' 35 
u - u '-0.32 


— 





ocr+0.54 
oo -0.93 


in q+67.8 

lz - -ll.l 


U.C. 


Abell 520 


14+0- 23 

u - ±4 -0.21 


PKS0745-191* 


19+ 010 
u - la -o.io 




1 


— 0.57 


8-0^ 5 5 ° 


> 0.89 


Abell 665 


O.SOlo 1 ;! 4 


IRAS 09104+4109 


21+ ' 26 
u - zi -0.24 







nq+0.36 
ao -0.80 


8.9+^ 


> 0.79 


Abell 773 


O.28I0.1 4 


Abell 963 


02+ ' 14 
u - uz -0.15 


— 





Q1+0.21 


— 


U.C. 


Abell 2163 


66+0 .13 
U - DD -0.12 


Zwicky 3146 


-0 19+ 011 
u - ly -o.io 


1 2+°' 3 
1- -0.3 





cq+0.17 
JO -0.26 


04 +i.6 


> 0.42 


Abell 2218 


n nc+0.14 
-O.O6To.x4 


Abell 1068 


85+ ' 15 
u - so -0.14 


q O + l.O 

o - z -0.8 


1 


1O+0.44 
iy -0.39 


1 7+ 11 


> 0.95 


Abell 2219 


59+0- 10 

u - oy -o.io 


Abell 1413 




4-1 





_i_n -^fi 
oo -0.19 


_I_Q a 

3.9+2 5 


> 0.89 


Abell 2319 


o.25iEJ:ii 


Abell 1689 


0.14+009 

— U.Uo 


4.1+" 

— l.o 





19+012 
ly -0.24 


16.2+i 2 5 6 


> 0.59 


AC 114 


o.5ii«: 22 


Abell 1704 


46+ ' 31 


q q + 5.2 

•"-1.6 





45 +i.38 


10.91JV 


> 0.79 






RXJ1347.5-1145 


n 49+O.16 
u -* y -0.15 


9 -+3.0 
z - ' -0.7 


1 


05+0- 17 


8-3l^ 3 2 


> 0.58 






Abell 1795 


05+ ' 05 


o 1+O.6 
z -!-0.4 





-17+0.11 
-I'-O.ll 


7 +1 -° 


> 0.76 






MS1358.4+6245 


S5+0' 32 
U.80_ 2g 


D -^-3.8 





QV+0.55 
y -0.22 


i c+1.9 
1 - 8 -0.7 


> 0.80 






Abell 1835 


n o 9 +0.09 
"•° -0.09 


q 9+16 
a - Z -0.8 





fi q+0.34 


3-i±!:S 


> 0.80 






MS1455.0+2232 


41+0 .17 
u - 41 -0.17 


4 1+ 4 - 3 

4 -1.6 





fi c+0.64 
DO -0.36 


i.8ti : ; 


> 0.63 






Abell 2029 


27+0' 04 
u - z '-0.05 


3 1+0- 5 
° -0.4 





45+0.10 
^ d -0.09 


6.61 2 ;! 


> 0.91 






Abell 2142 


S5+ ' 10 


3 6+ 18 
°- D -i.o 





qc+0.27 
OJ -0.08 


15-0lli 5 4 


> 0.80 






Abell 2204 


0.611°;- 


q q+0.9 

°-°-0.8 


1 


9+0.37 
UZ -0.38 


1 9+ 14 

-0.7 


> 0.93 






Abell 2261 


0.84±°;1? 


3.811-1 


1 


r-i+0.39 
ol -0.61 


5.9 +3 ' 7 


> 0.84 






MS2137.3-2353 


59+0- 29 
U.O»_ 28 







60+0- 39 
DU -0.19 


3.3l 4 ^ 


> 0.87 






Abell 2390 


43+0- 32 


2.9+™ 


1 


17 +0.42 

^-'-o.ss 


5 5+14.9 


> 0.57 






MEAN 


0.35 ±0.30 


3.4 ±1.3 


0.72 ±0.39 


6.3 ±4.7 




MEAN 


0.44 ± 0.39 



A summary of the results on intrinsic absorption from the ASCA data. Columns 2 — 5 list the differences between the spectrally- 
determined column densities and the nominal Galactic values for spectral models B, C, D and C Also listed are the (90 per cent 
confidence) constraints on the covering fraction of the intrinsic absorber, /, using spectral model C. All results assume solar mctallicity 
in the absorbing material. For PKS0745-191* , the ASCA data indicate a Galactic column density 3.5 X 10 21 atom cm -2 which is lower 
than the nominal value of 4.24 X 10 21 atom cm~ 2 (Dickey & Lockman 1990). The lower value of 3.5 X 10 21 atom cm~ 2 has been used in 
the calculation of the intrinsic column densities. The errors on the mean values are the standard deviations of the distributions. 
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Table 10. The masses of intrinsic X-ray absorbin 



g gas 



Cluster 



^abs 

(kpc 



M abs [ZQ 
(10 n M e 



A^abs [-^ICm] 



M(r < r abs ) 
(10 9 M yr" 1 ) 



■^acc 

(10 n M r . 



Abcll 478 
PKS0745-191 
IRAS 09104+4109 
Zwicky 3146 
Abell 1068 
Abell 1689 
Abcll 1704 

RXJ1347.5-1145 
Abcll 1795 

MS1358.4+6245 
Abcll 1835 

MS1455.0+2232 
Abcll 2029 
Abell 2142 
Abell 2204 
Abcll 2261 

MS2137.3-2353 
Abcll 2390 



120^ 

iiftg 

i qc+38 
1JO -39 
123^ 

oo+69 
»»_19 

114+ 90 

1 24 +22 
lz -29 

79+65 

' y -50 

1261 17 

qo+32 

1361™ 
94l? 6 

62±~ 



7.30] 
13.31"' 



5.52] 
18.93^ 



-,+4.60 
-2.61 
1+12.72 
-8.30 
14 02+ 27,83 

,+4.83 
-1.94 
1-13.13 
-9.60 
on 1 q+22.05 

^u.id_ 11 g9 
10.50l 24 9 5° 
11.39™ 

12 g7 +35.97 
lz - y '-11.32 
17 29+ 17 ' 27 

43.60+ 27 ' 57 

15.98 . 

10-Hl 6 .92 
19.8ll« 7 26 



30.81 
+4.35 
-4.17 
-9.06 



10.90 

58.16 
3.62 



+ 10.27 
10.51 
+66.00 
-38.17 
11.28 
2.40 



14.60l 9 - 2 2 

26-61™ 
28.051^-66 

11 no+9.67 
ll.UJ_3 gg 

37.86l 2 S' 2 8 
40.27l«; 9 ° 
21.00lf 9 °° 
22.78l^; 2 { 
20.96±|;|| 
25 93+ 71 - 94 
34.57l^;« 
87.201^-14 

2o-2i±ii:J| 

39 fi'?+36.51 
21.80l«°;g* 

lie OO + 132.00 
llD.OZ_ 76 34 

7 24 +22.57 
'• Z4 _4.79 



6281™ 
909l 344 
9111 292 
1050lf 4 9 

7fi o+232 

g62 +890 
5DZ _347 

1440l 7 ° 8 
645l» 7 

232I133 
1430™ 
15961™ 
839+ 152 
374+ 219 

13371^5 
329 +20l 
•J zy _301 

293l^ 3 



15-70±i; 7 ° 
22.731^0 
22.77+JI 

26.25 ^ 

19.231^° 

91 kc;+22.25 

zi.oo_ 8 6g 
9.65™ 
36-OOtl™ 
16 13+ 3 - 43 

5.80l 20 3 27 

qtr 7 c-+32.67 

39.90l£ 73 
20.98l 3 ;«° 

33.43™ 
8.23l«° 3 
41.08™ 

7 qo+14.33 

<-J-5_2.47 



Notes: r a b s values are the radii at which the cooling time of the gas first exceeds 5 X 10 9 yr, determined from the absorption-corrected 
deprojection results. M a b s values are the masses of absorbing material within radii r a b s implied by the observed intrinsic column densities 
(determined with spectral model C) assuming the absorption to be due to cold gas with solar metallicity (column 3) or Z = 0.4Zq 
(column 4). Column 5 lists the integrated mass deposition rates (corrected for the effects of absorption) within radii r a b s . M acc values 
are the masses of material accumulated within radii r a b s during the first 5 Gyr. Errors incorporate the statistical uncertainties in both 
the intrinsic column densities (determined with spectral model C) and the deprojection results. Errors on the r a b s values are the 10 and 
90 percentile values from 100 Monte Carlo simulations. Errors on the M a b s values are the extrema incorporating the joint errors in r ab s 
and AJV a b s . Results are not tabulated for Abell 1413 since the central cooling time in this cluster exceeds 5 Gyr. 



© 0000 RAS, MNRAS 000, 000-000 



The properties of cooling flows in X-ray luminous clusters of galaxies 35 

Table 11. The reprocessed luminosities and observed IRAS fluxes 



ircpro. S%0 S\00 AR £l-1000 M m 

Cluster (lO^ergs- 1 ) (mjy) (mjy) (arcmin) (lO^crgs- 1 ) 



Abell 478 


10.5 


570 ± 58 


5490 ± 440 


0.55 


60 ± 19 


PKS0745-191 


13.4 


< 294 


< 1680 




< 29.0 


IRAS 09104+4109 


13.3 


560 ± 39 


< 324 


— 


< 390 


Zwicky 3146 


6.15 


< 111 


880 ± 56 


1.53 


< 111 


Abell 1068 


9.05 


680 ± 38 


950 ± 108 


0.02 


59 ± 18 


Abell 1413 


5.71 


< 111 


< 333 




< 14.2 


Abell 1689 


3.46 


< 141 


< 357 




< 27.3 


Abell 1704 


6.84 


90 ±41 


280 ± 83 


0.43 


27 ± 11 


RXJ1347.5-1145 


49.7 


< 141 


< 498 




< 196 


Abell 1795 


1.95 


< 144 


< 354 




< 3.23 


MS1358.4+6245 


6.45 


< 81 


< 222 




< 51.6 


Abell 1835 


16.7 


140 ± 64 


< 504 


0.05 


< 61.5 


MS1455.0+2232 


12.3 


< 96 


< 378 




< 46.6 


Abell 2029 


5.13 


< 96 


< 351 




< 3.98 


Abell 2142 


5.47 


130 ±51 


480 ± 156 


2.50 


7.2 ± 2.9 


Abell 2204 


19.2 


< 288 


1760 ± 306 


2.58 


< 64.8 


Abell 2261 


14.9 


150 ± 30 


410 ± 53 


1.66 


45 ± 15 


MS2137.3-2353 


11.0 


< 69 


< 249 




< 46.9 


Abell 2390 


17.4 


130 ± 60 


750 ± 207 


1.73 


66 ± 26 



Notes: Column 2 lists the X-ray luminosities expected to be reprocessed into far infrared emission by the intrinsic absorbing material 
in the clusters. The best-fit parameters determined with spectral model C have been used in the calculations. Columns 3 and 4 lists 
the observed 60 and lOOfim IRAS fluxes within a 4 arcmin (radius) aperture centered on the X-ray centres (Table 3), determined with 
the IPAC SCANPI software and co-added IRAS scans (we use the median of the coadded scans). Error bars are the root- mean-square 
deviations in the residuals, external to the source extraction regions, after baseline subtraction. Where no clear detection was made, an 
upper limit equal to three times the r.m.s. deviation in the residuals is given (although systematic errors could exceed these limits in 
individual cases). Where a possible detection is made, Column 6 lists the in-scan separation, AR (in arcmin), between the peak of the 
100/im emission (or 60^m emission in the case of Abell 1835) and the X-ray centre. Column 7 lists the total 1 — 1000/im luminosities 
calculated using equation 4. 
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